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1.  General ■  The  subject  grant  supported  research  work  in  control  of 

distributed  parameter  systems,  control  of  nonlinear  systems,  and  the 
mathematical  modelling  of  such  systems  by  the  Principal  Investigator 
and  his  research  assistants  during  the  period  Feb.  15,  1934  -  Oct. 
15,  1985.  In  addition  to  salary  support  for  the  Principal  Investi¬ 
gator  and  his  research  assistants,  funds  were  used  to  bring  visitors 
to  the  University  of  Wisconsin  campus  as  scientific  consultants  for 
short  periods,  to  support  scientific  computing  relevant  to  the  re¬ 
search  program,  to  support  domestic  and  foreign  travel  by  the  Prin¬ 
cipal  Investigator  and  to  purchase  needed  equipment  and  supplies  for 
the  UW  MIPAC  (Modelling,  Information  Processing  and  Control)  Facil¬ 
ity,  for  which  the  Principal  Investigator  serves  as  coordinator. 
Below  we  describe  these  activities  in  greater  detail.  A  copy  of  a 
recent  research  article  is  attached  as  an  Appendix  to  this  report. 
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2.  Areas  of  Research  Emphasis.  During  the  grant  period  a  wide  rang¬ 
ing  program  of  scientific  research  in  mathematical  systems  theory 
was  carried  forward.  In  the  rest  of  this  section  we  summarize  prog¬ 
ress  in: 


(a)  Control  Theory  of  Distributed  Parameter  Systems; 

(b)  Control  of  Nonlinear  Systems; 

(c)  Modelling  of  Distributed  Parameter  Systems; 

(d)  Coefficient  Identification  in  Distributed  Systems. 

(a)  The  Principal  Investigator  has  worked  and  supervised  research 
in  the  area  of  control  of  distributed  parameter  systems  for  many 
years,  developing  theories  applicable  to  wave  (hyperbolic)  processes 
and  diffusion  (parabolic)  processes  as  well  as  many  other  areas  im¬ 
portant  in  applications.  During  the  grant  period  two  research  assi¬ 
stants  were  supervised  in  carrying  out  doctoral  dissertation  work  in 
these  areas. 

The  first  of  these  assistants,  Richard  Rebarber,  received  the 
PhD.  degree  in  August,  1984.  His  dissertation  was  concerned  with 
control  canonical  forms  and  spectral  assignment  problems  related  to 
control  of  infinite  dimensional  systems  generically  described  by 


x  =  A  x  +  B  u 


wherein  the  operator  A  is  the  generator  of  a  holomorphic  semigroup. 
Because  of  recent  experimental  evidence  and  related  mathematical 
modelling  work,  which  we  describe  under  (iii)  below,  it  is  now  clear 
that  the  most  realistic  models  for  vibration  of  beams,  plates,  etc., 
are  actually  systems  of  this  sort,  involving  frequency-proportional 
damping.  Rebarber's  results  deal  with  the  extent  to  which  the  closed 
loop  spectrum  of  the  system  obtained  with  linear  feedback  u  =  K  x 
can  be  specified  by  choice  of  the  feedback  operator  A. 
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A  second  assistant,  Katherine  Kime,  is  expected  to  receive  the 
PhD.  degree  in  August,  1986.  During  the  grant  period  she  was  parti¬ 
ally  supported  with  grant  funds  in  her  research  on  control  of  elect¬ 
romagnetic  fields,  satisfying  the  Maxwell  equations,  by  means  of 
electric  currents  flowing  on  the  boundary  of  the  spatial  region  in 
which  the  fields  are  defined.  This  research  has  extended  earlier 
work  carried  out  by  the  Principal  Investigator  under  AFOSR  auspices. 
Ms.  Kime  has  been  able  to  demonstrate  the  controllability  of  finite 
energy  states  in  a  three  dimensional  spherical  region  and,  with  add¬ 
itional  smoothness  requirements,  in  more  general  regions.  We  ex¬ 
pect  the  work  ultimately  to  have  implications  for  design  of  radar 
non-reflecting  surfaces,  etc. 


(b)  We  have  been  interested  for  some  time  in  the  control  of  non¬ 
linear  systems  exhibiting  self-excited  oscillations  because  of  their 
importance  in  the  study  of  flutter  phenomena  in  aircraft  design. 
Emphasis  has  been  placed  on  systems  of  the  form 

x  =  Ax  +  Bu  +  Cy 
y  =  g  ( x ,  u ,  y ) 

wherein  x  is  the  state  of  an  elastic  system  with  basically  linear 
dynamics  and  y  is  the  state  of  a  (generally  lower  dimensional)  non¬ 
linear  system  interracting  with  the  elastic  system  to  produce  non¬ 
linear  oscillations  whose  amplitudes  may  become  large  as  their  fre¬ 
quency  approaches  one  of  the  resonant  frequencies  of  the  elastic 
structure.  We  are  concerned  not  only  with  the  question  of  using  the 
control  input,  u,  to  alleviate  the  effects  of  the  self-excited  osci¬ 
llations  in  the  presence  of  complete  information  about  the  total 
system  state  (x,y),  but  also  with  the  question  of  system  state  esti¬ 
mation  and  resultant  control  specification  when  a  lower  dimensional 
output 


u  =  d(x,y) 


is  all  that  is  available.  A  doctoral  candidate,  Thomas  Svobodny, 
was  supported  during  the  grant  period  in  the  initial  stages  of  his 
work  in  connection  with  this  latter  question.  He  has  been  seeking 
to  extend  to  nonlinear  systems  exhibiting  self-excited  oscillations 
some  work  done  earlier  by  the  principal  investigator  concerning 


adaptive  rejection  of  periodic  disturbances,  v(t),  in  a  system 

x  =  Ax  +  Bu  +  v 

Promising  initial  results  on  observability  of  linear  systems  with 
time-periodic  coefficients  are  being  extended  to  be  applicable  to 
the  nonlinear  systems  of  interest.  We  have  also  been  acquiring  lab¬ 
oratory  data  on  actual  oscillating  systems  to  serve  as  test  cases 
for  the  identification,  estimationand  control  procedures  being  deve¬ 
loped  in  these  researches. 

(iii)  During  the  contract  period  our  main  mathematical  modelling  re¬ 
search  activity  has  been  in  the  area  of  development  of  mathematical 
models  to  replicate  the  observed  f requency-proportional  "structural" 
damping  properties  of  actual  beams.  While  this  work  has  been  comp¬ 
leted  under  a  successor  grant,  AFOSR  85  -  0283,  we  can  state  unequi¬ 
vocally  that  the  major  conceptual  advances  were  made  under  the  aegis 
of  the  subject  grant,  AFOSR  84  -  0088,  during  the  period  when  the 
Principal  Investigator  was  visiting  the  University  of  Florida  at 
Gainesville.  Following  many  unsuccessful  attempts,  a  model  in  the 
form  of  an  integro-partial  differential  equation 

P3—  +  2y  f  h(x,  £ )  fSfn_(X,  t)  -  2fw_($,t)]d$  +  EI^J  =  0 
9t^  J0  ^tSx  Bt'Bx  J  3x* 

was  developed  which,  with  y  and  the  "interraction  kernel"  h(x,$) 
correctly  chosen  relative  to  the  mass  density  p  and  bending  modulus 
El,  results  in  very  close  spectral  matching  with  experimental  data 
taken  from  various  types  of  beams  excited  into  vibratory  motion  in 
a  laboratory  setting.  We  have  high  hopes  that  this  model  will  find 
wide  application  in  theoretical  damping  considerations  relative  to 
large  space  structures  and  other  similar  topics  of  current  interest. 

( iv)  We  have  also  been  interested  in  the  problem  of  identifying  the 
coefficients  of  a  partial  differential  equation,  assumed  to  lie  in  a 
particular  class  of  such  equations,  from  data  taken  from  solutions 
of  the  equation.  Such  problems  arise  in  porous  flow  problems,  where 
the  potential  function  <P  can  be  measured  in  the  field  and  is  assumed 
to  satisfy  an  elliptic  equation  of  the  form 
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V'  (P7<?)  =  f 

with  p  the  permeability  of  the  flow  medium  and  f  the  source  distri¬ 
bution  function.  A  variety  of  L1  and  L”°  methods  admitting  approxi¬ 
mate  realization  in  linear  programming  codes  have  been  studied  and 
tested  computationally  with  very  promising  results.  A  graduate 
research  assistant,  Robert  Acar ,  has  been  working  with  the  Princi¬ 
pal  Investigator  in  this  area. 

A  second  type  of  identification  problem  concerns  equations  of 
evolution  type  describing  vibrations  of  physical  continua,  such  as 
the  wave  equation 


r32w  _  9  f  k9w  -j 
9t2  3x*-  3x  ' 


=  0 


or  the  beam  equation  (undamped  here) 


Er- 


9x 


0  . 


It  is  assumed  that  a  scalar  functional  of  the  state  is  recorded, 
whose  power  spectrum  reveals  the  natural  frequencies  of  vibration 
of  the  system.  The  problem  is  to  recover  whatever  information  is 
recoverable  about  the  coefficient  functions  of  the  system  from  the 
spectral  data  at  hand.  A  perturbation  procedure  based  on  spanning 
and  independence  properties  of  the  squared  eigenfunctions  has  been 
initially  formulated  and  will  be  a  subject  for  future  study. 


\ 
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3 •  Other  Scientific  Activities  under  AFOSR  84  -  0088 . 

(a)  Visitors_. _  Grant  funds  were  used  to  bring  a  number  of  visitors 

to  the  U.W.  campus  for  short-term  scientific  consulting  purposes. 
These  visitors,  with  the  approximate  duration  of  their  stay,  their 
home  affiliation  and  their  area  of  expertise,  included: 

E.  Fernandez  (3  weeks)  VP I,  Functional  Equations; 

K.  Kunisch  (3  days)  Univ.  of  Graz,  Identification; 

Q.  Iqbal  (3  weeks)  Univ.  of  Karachi,  Algebraic  Systems. 

(b)  Travel^  Grant  funds  were  used  to  support  domestic  travel  by  the 

Principal  Investigator  to  the  1984  Control  and  Decision  Conference 
in  San  Antonio,  Texas  and  to  AFOSR  headquarters  in  the  spring  of 
1985.  Funds  also  supported  attendance  by  the  Principal  Investigator 
and  his  assistant,  Robert  Acar ,  at  the  October,  1985  meeting  on 
Systems  Identification  and  related  topics  at  the  University  of  Okla¬ 
homa  in  Norman,  OK.  Foreign  travel  sponsored  by  the  grant  consisted 
of  the  July,  1984  trip  to  Vorau  and  Graz,  Austria,  for  the  1984 
Vorau  Conference  on  Identification  and  Control  of  Distributed  Para¬ 
meter  Systems  and  for  a  week  of  scientific  meetings  in  Graz  with 
Professors  Kappel  and  Kunisch. 

(c)  Scientif ic_Computingi  Supplies^  Equipment.  Grant  funds  were  used 

to  purchase  time  for  scientific  computing  on  the  UW  MACC  Univac  1110 
Computer  for  use  by  the  Principal  Investigator  and  his  assistants. 
Funds  were  also  used  for  various  categories  of  supplies  and  office 
needs  such  as  stationery,  postage,  long  distance  telephone  calls, 
etc.  A  variety  of  equipment  to  be  used  in  the  construction  of  labo¬ 
ratory  models  for  study  and  comparison  with  mathematical  models  was 
purchased  and  put  into  use  in  the  UW  MIPAC  Model  Development  Unit  at 
1307  University  Ave . ,  Madison.  Certain  modifications  to  MIPAC's 
HP5451C  Fourier  Analyzer  were  paid  for  with  grant  funds  and  a  vari¬ 
ety  of  related  supplies  such  as  computer  discs  and  tapes,  graphics 
supplies,  model  construction  supplies,  etc.  were  paid  for  with  grant 
funds.  Details  will  be  provided  in  the  Financial  Report  prepared 
by  our  Office  of  Research  Administration,  Financial. 


Grant  funds  sponsored  research  resulting  in  a  number  of  scien¬ 
tific  reports  issued  over  the  lifetime  of  the  grant  and  paid  for 
their  preparation  and,  in  some  cases,  for  journal  reprints.  A  list 
of  these  publications  follows: 

1.  On  the  Dir ichlet-Neumann  boundary  control  problem  associated 
with  Maxwell's  equations  in  a  cylindrical  region.  To  appear 
in  SIAM  Journal  on  Control  and  Optimization,  1986. 

2.  Frequency/period  estimation  and  adaptive  rejection  of  periodic 
disturbances.  To  appear  in  SIAM  Journal  on  Control  and  Optimi¬ 
zation,  1986. 

3.  A  Floquet  decomposition  for  Volterra  equations  with  periodic 
kernel  and  a  transform  approach  to  linear  recursion  equations. 
Submitted  to  the  Journal  of  Differential  Equations. 

4.  Mathematical  models  for  the  elastic  beam  and  their  control- theo¬ 
retic  implications .  Proc.  1984  Autumn  College  on  Semigroups  and 
their  Applications,  Int'l.  Centre  for  Theoretical  Physics, 
Trieste,  Italy,  November  1984.  To  be  published  by  Birkhauser. 

The  above  are  by  the  Principal  Investigator.  Also  sponsored  by  the 
grant  was  the  following,  by  Richard  Rebarber: 

5.  Control  canonical  forms  and  spectral  assignment  for  holomorphic 
semigroups.  Thesis,  University  of  Wisconsin,  Madison,  Aug.  1984 
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APPENDIX : 

"Mathematical  Models  for  the  Elastic  Beam 
with  Frequency-Proportional  Damping" 


by 

David  L.  Russell 


University  of  Wisconsin,  Madison 


Introduction  to  Elementary  Theories  of  Elastic  Beam  Motion 


The  less- than-dashing ,  rather  pedestrian  term  "beam"  has  a 
guaranteed  soporific  effect  on  all  but  the  more  dedicated  devotees 
of  the  science  of  structural  mechanics.  But  it  refers,  of  course, 
to  the  indispensable  linear  structural  elements  (those  seeking  new 

terminology,  please  take  note!)  without  which  many  complex  constr¬ 
uctions  would  be  infeasible.  Consequently,  it  is  not  surprizing 
that  mathematical!  models  for  the  elastic  beam  have  an  ancient  and 
honorable  pedigree  attested  to  by  the  attentions  of  some  of  the 
most  honored  historical  and  current  figures  of  applied  mathematics. 

The  Euler  -  Bernoulli  equation  for  the  motion  of  thin  elastic 

beams  is  at  least  two  hundred  years  old.  Denoting  the  mass  density, 
per  unit  length,  by  p(x)  and  the  second  moment  of  the  modulus  of 
elasticity  about  the  elastic  axis  (about  which  the  first  moment  of 
the  modulus  of  elasticity  vanishes)  by  EI(x),  0  £  x  £  L,  it  is  as¬ 
sumed  that  the  energy  associated  with  motion  in  the  x,w  plane,  in 
which  the  elastic  axis  is  given  by  w  =  0,  can  be  adequately  repre¬ 
sented  by 


£(w. 


3w .  _ 


r  f  3w 

l  i,  J 


+  El 


n 


If  it  is  assumed  that  no  work  is  done  on  the  beam,  either  by  exter¬ 
nally  applied  forces  or  by  internal  dissipative  mechanisms,  an  easi 
application  of  the  principle  of  virtual  work  shows  that  the  motion 
must  be  governed  by  the  partial  differential  equation 


.  92w  a2 

H - -5  + - 

at2  ?x 


a  f 

^  l 


^  ]  =  0 
ax-  J 


(1.2) 


together  with  "natural  boundary  conditions"  (four  are  required  in 
all)  implying  that 


-w  alw_  j L 
: 2  at?x  1 0 
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A  variety  of  conditions,  corresponding  to  different  kinem.atic  and 
dynamic  assumptions,  suffice  for  this  purpose.  The  "cantilever" 


configuration,  wherein  the  beam  is  clamped  at  one  end,  say 
and  left  free  at  the  other,  x  =  L,  is  described  by 


0, 


w(0,t)  =  0,  2^(0, t)  =  0, 

Px 


(i 


"W 


2 


(  L  ,  t )  =  0  ,  (El  &S) 

9>:  9:X 


=  0 


x  =  L 


If  a  lateral  force,  V,  and  moment,  y ,  act  at  x  =  L,  then  the  righ 
hand  sides  of  the  equations  in  (1.4)  are  replaced  by  and  9,  res 
pertiveiy.  We  will  describe  some  other  cases  subsequently. 


The  F.ayleigh  beam  model  is  a  fairly  modest  modification  of  t 


Euler  -  Bernoulli  model,  dating  back  some  hundred  years  or  so.  T 
obtain  it  the  beam  element  of  width  centered  at  x  is  additional 


endowed  with  mass  moment  of  inertia 


i  I  (x)  .  The  term  I 


is  added  to  the  energy  integrand  of 
partial  differential  equation 


(1.1)  and  there  results  the 


3.  f  i  _3f«_ 

9;:  ^  p  9t29x 


1  . 


with  corresponding  modifications  of  the  boundary  conditions  which 
we  will  not  detail  at  this  point. 


originally  rectangular  w,x  cross  section  is  transformed  into  a  par 
allellogram  as  shown  in  Figure  1(b).  ~.ie  total  rotation  angle  of 
an  element  is  thus 

*  ( x ,  t )  -  6  ( ,  t )  +  f  ( >: ,  t ) 

where  for  small  motions  the  bending  angle  9  is  given  approximately 
by 

c(x.t)  :  (x.t) 

9x 

and  £  is  the  angle  of  shear.  Consequently 

p(x,t)  2  t(x,t)  -  (x,t)  . 

3x 


The  energy  expression  used  is 


C  ( w ,  -  -  ,  * ,  -  ' 
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ax 
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where  K(x)  is  the  modulus  of  elasticity  in  shear.  Again,  application 
of  the  principle  of  virtual  work  quickly  leads  to  the  Timoshenko 
equations 


a2w 


P  -  2 
at2 


^  f  t,- 

'■  t  K  ' 


ax1 


ax 
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'  t  _  aw  •) 
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ax 
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As  inspection  quickly  shows,  these  equations  have  the  form  of 


coupled  wave  equations.  In  general  the  wave  speeds  f  K  ( x )  /  p  ( x 

1  1/2 

enn  |  r.  1  ( x ) I  ^  ( x )  j  are  different,  resulting  in  a  hyperbolic 

tern  with  four  families  of  characteristic  curves,  consisting  of 
pairs,  each  pair  corresponding  two  waves  of  the  same  speed  mov 
opposite  directions.  The  boundary  conditions  come  out  of  the 
calculations  as  produce  the  equations  (1.6)  and  (1.7).  In  the 
case  of  the  cantilever  configuration  these  conditions  are  easi 
seer,  to  take  the  form 


two 

i  / ; 
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sys- 

two 

ing  in 


same 


-  i 


K ( L )  j  ( L , t 

)  -  +(L,t)  ] 

1  ax 

J 

where,  again. 

and  p  are 

respect ively , 

at  x  =  L . 

In  the  constant  coef 

w  (  0  ,  t )  =  0,  H-1  (  0  ,  t )  =  o, 

<p(t),  E I  ( I )  —  (  L  , 

ax 


1.8) 

1.9) 


comparison  with  K,  as  tends  to  be  the  case  for  thin  beams,  the  first 
equation  can  be  differentiated  with  respect  to  x,  solved  f 
substituted  into  the  second  equation  to  obtain 


a2*  . 

or  ana 

9x  2 


2 .  Internal  Damping  Mechanisms . 

In  traditional  applications  it  is  possible  to  ignore  internal 
damping  effects  to  a  large  degree.  This  is  the  case  because,  in 
those  applications  interest  is  largely  confined  to  static  deforma¬ 
tion  or,  at  a  slightly  more  sophisticated  level,  to  the  vibrations 
produced  by  periodically  imposed  loads,  such  as  might  be  produced, 
for  example,  by  a  train  of  heavy  cars  crossing  a  bridge  or  a  power¬ 
ful  rotating  engine  mounted  on  an  aircraft  wing.  Damping  effects 
in  these  cases  can  be  largely  ignored  because  the  standard  rule  of 
design  is  to  keep  the  natural  frequencies  of  vibration  of  the  struc¬ 
ture  well  separated  from  the  frequencies  of  expected  periodic  loads. 

As  aircraft  have  grown  ever  more  sophisticated  and,  particular¬ 
ly,  as  the  construction  of  platforms  in  space  has  been  contemplated , 
it  has  become  increasingly  necessary  to  take  internal  energy  losses 
into  account.  As  structures  have  grown  lighter  and  stronger  they 
have  also  become  more  elastic  and  hence  more  subject  to  vibration 
problems.  At  the  same  time,  pointing  accuracy  and  platform  stabil¬ 
ity  requirements  have  made  the  suppression  of  such  vibrations  abso¬ 
lutely  essential.  It  has  become  necessary  to  augment  natural  damp¬ 
ing  by  the  introduction  of  active  control  mechanisms  and/or  care¬ 
fully  matched  passive  damping  systems.  Paradoxically,  as  more  elas¬ 
tic  structures  with  less  inherent  damping  have  evolved,  and  as  ac¬ 
tive  control  has  increasingly  been  brought  to  bear,  the  importance 
of  a  thorough  understanding  of  natural  damping  mechanisms  and  of  a 
capability  to  accurately  estimate  its  effects  has  become  more, 
rather  than  less,  essential.  This  is  true  because  artificial  damp¬ 
ing  achieved  through  active  control  implementation  is  inevitably 
limited  in  band  width;  at  higher  frequencies  its  effectiveness  de¬ 
creases  or  disappears  entirely  -  indeed,  phase  lags  may  result  in 
positive,  rather  than  negative,  work  being  done  by  the  control 
forces  on  the  higher  frequency  modes.  That  the  controls  used  are 
generally  effective  is  largely  thanks  to  the  fact  that  natural  damp¬ 
ing  mechanisms  tend  to  become  stronger  with  increasing  frequency, 
so  that  artificial  controls  can  concentrate  on  suppression  of  the 
modes  associated  with  lower  frequencies,  leaving  higher  frequency 
"spillover"  effects  to  be  "mopped  up"  by  the  strong  internal  damping 
operative  at  the  higher  end  of  the  spectrum.  As  internal  damping 
decreases,  it  is  ever  more  important  to  be  able  to  estimate  its 
precise  strength  so  that  the  band  width  requirements  for  the  applied 


active  controls  can  be  appropriately  specified. 


That  damping 
been  understood^, 


rates  tend  to  increase  with  frequency  has  long 
Thus  simplistic  models  for  damping  such  as 


■s 

0 


0  , 


arc  yererolly  w  • -L  *-  y  «.».■/«. 

which  produce  uniform  damping  rates,  c-annot — fee — seri^ncl-y 
The  recognition  that  damping  rates  in  beams  increase  with  increasing 
frequency  goes  back  at  least  to  Lord  Kelvin  in  Britain  and  Robert 
Voigt,  a  distinguished  German  physicist,  both  working  at  the  end  of 
the  last  century  [  fs|  ].  The  Kelvin  -  Voigt  damping  model  applies,  i 
principle,  to  the  vibrations  of  any  linear  elastic  system.  The 
model  hypothesizes  that,  whatever  the  linear  operator  describing 
the  elastic  restoring  forces  may  be,  the  damping  forces  involve  a 
positive  multiple  of  that  operator  acting  on  the  system  velocity 
rather  than  the  displacement.  This  theory,  modifying  the  basic 
Euler  -  Bernoulli  model,  yields  an  eqxiation  which  may  be  taken  to 
have  the  form 


n‘w 
P  2 
9t 2 


+  2yP  ~—~2 

f  EI  ^2  1 

f  EI  1 

=  0 

at9x2 

L  ax2  J 

9x2 

1  ax2  j 

with  appropriately  modified  boundary  conditions.  If  the  fourth 
order  elasticity  operator  is  denoted  by 


Aw 


i  a' 


ax' 


EI 


a^w  i 

ax2  J 


then  the  operational  form  of  the  equation  is 


+  2y  A—  +  Aw  =  0 
dt 2  dt 


If  the  eigenvalues  of  the  positive  self-adjoint  operator  A  are  , 

so  that  the  natural  frequencies  of  the  undamped  system  are  uj.  = 

i/2 

(X^)  .  then  the  damped  system  (2.1)  may  be  seen  to  have  exponen¬ 

ts. .t 

tial  solutions  e  sfv  ,  where  <p,  is  the  corresponding  eigenvector  of 


A  and  o-  satisifes  the  quadratic  equation 


-k2  +  2^wk2®k  +  uk2  =  0  • 


Thus  c-  is  given  by 


Cu  —  >u 


Jk  1  l  r  °-k 


(2.2) 


Assuming  >  to  be  small,  we  see  that  the  Cj.  are  complex  for  some 

finite  number  of  values  of  k  with  the  damping  rate  proportional 
to  the  square  of  the  frequency.  Critical  damping  occurs  for  u  = 
l/V;  for  uj_  larger  than  this  the  modes  are  overdamped  with  one  of 

the  values  given  by  (2.2)  going  to  -*•  and  the  other  tending  to 

the  value  The  values  Cj.  ,  in  fact,  lie  in  the  locus  shown 

in  Figure  2.1,  consisting  of  a  circle  of  radius  — -  ,  centered  at 

the  point  ( ,  0 )  and  the  portion  of  the  negative  real  axis  to  the 

left  of  (~",0).  For  very  small  values  of  y  all  we  would  expect 

to  see  would  be  the  quadratic  dependence  of  the  damping  rate  on  the 
frequency.  Whether  the  overdamping  predicted  by  the  model  has  ever 
been  observed  in  the  laboratory  is  unknown  to  this  writer  but  it 
seems,  on  the  face  cf  it,  to  be  unlikely. 


•  %  ^  N  '  "L'  »  '  *  ’  V'  *L  «  '  \  *  ■  *  »*•  •  '  ,  4.’  .  '  v"«  ’  •  *  -  '  *'  '  ^  '  •.* 
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Fig.  2.1 


In  any  case,  as  we  will  document  in  the  final  section  of  this 
paper,  our  experiments  at  UW  MIPAC,  in  agreement  with  experimental 
data  acquired  by  many  other  researchers,  indicates  a  predominantly 
linear  relationship  between  the  damping  rate  and  the  frequency;  a 
behavior  widely  known  as  structural  damping..  We  will  also  see,  both 

from  the  experimental  and  theoretical  point  of  view,  that  damping 
of  this  sort  is  by  no  means  a  "rule  of  nature";  indeed  there  is 
some  reasons  to  believe  that  it  may  only  apply  to  thin  beams  for 
which  the  Euler  -  Bernoulli  model  is  adequate  in  the  conservative 
case.  More  on  this  later. 

There  is  no  difficulty  in  obtaining  mathematical  models  which 
exhibit  this  sort  of  damping/frequency  relationship.  Representing 
an  arbitrary  linear  oscillator,  without  damping,  by 


dfw 

dt2 


+ 


Aw 


0  , 


the  simplest  mathematically  viable  example  of  a  system  exhibiting 
structural  damping  behavior,  treated  extensively  by  the  author  and 
G.  Chen  in  [  E  ] ,  is  given  by 
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I 

— ~  +  2y  Al/2w  +  Aw  =  0  ,  (2.3 

;  dt2 

where  A17*2  denotes  the  positive  square  root  of  the  positive  self- 
-ad joint  operator  A  and  y  >  0.  Attempting  a  solution  of  the  form 
okt 

e  <Pk  again,  we  find  now  that 


ck  =(->■  +  {y2  -  l) 1/2 )  uk 

so  that,  for  0  <  y  <  1  the  ok  lie  on  the  rays  in  the  left  half 
plane,  shown  in  Figure  2.2,  making  an  angle 


a  =  tan  1  [y/  (1  -  v2)] 


To  explain  why  we  cannot  embrace  the  model  (2.3)  with  unre¬ 
strained  enthusiasm,  let  us  consider  the  constant  coefficient  case 
wherein 


%  '■>  . 
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Aw  =  ei  alw 

p  3x4 


(2.4; 


Let  us  denote  boundary  conditions  for  this  operator  generically  by 


<  &  ,  w  >  =  0 , 


(2.5 


where  £  is  a  linear  functional  defined  on  the  domain  of  A.  We  will 
say  that  the  boundary  conditions  given  for  A  are  paired  if  for  each 

boundary  condition  (2.5)  appearing  in  the  set,  the  condition 


<  £  ,  w' '  >  =  0 


(2.6) 


appears  in  the  set  also.  Examples  include  hinged  endpoints,  for 
which  (taking  the  endpoint  in  question  to  be  x  =  C  ) 


w  (  0 ,  t ) 


=  0,  EI  (0,t)  =  0, 

9::2 


and  freely  clamped  endpoints  for  which  (constant  coefficient  case; 

(  0  ,  t )  =  C  ,  -  EI  (  0  ,  t )  =  0  . 

3x  9xJ 


For  paired  boundary  conditions  the  positive  square  root  of  the 


a4w 


a2w 


fourth  order  operator  — ,  is  -  2__  and,  redefining  y  slightly,  the 

3x*  9x^ 

modified  Euler  equation  in  this  case  takes  the  form 


32w  „  g3w 

p  —~2  -  2y 

St'2  dtdx2 


+  EI  =  0  . 

ax'* 


2.7) 


In  this  equation  the  damping  term  is  rather  easy  to  understand  from 
the  physical  point  of  view  -  it  is  a  lateral  force  acting  on  the 
beam  which  is  negatively  proportional  to  the  bending  rate. 


In  those  cases  where  the  boundary  conditions  for  A  are  not 
paired,  as  in  the  case  of  clamped  or  free  endpoints,  respectively 
corresponding  to 


w ( 0  ,  t )  =  0, 


(0,t) 


0, 


cr  w 

9>: 


El  2.*”  (  0  ,  t )  =  0, 


El 


(0,t) 


0, 


for  example,  the  positive  square  root  of  the  fourth  order  operator 
is  not  the  negative  second  order  operator  and  the  nature  of  the 
damping  term  in  (2.3)  is  such  as  not  to  admit  a  ready  physical  in- 
terpretation  in  terms  of  the  properties  of  standing  waves  in  the 
beam.  If  the  motion  is  conceived  of  in  terms  of  travelling  waves 
reflecting  from  the  boundaries  of  the  beam  then,  since  the  speed  of 
travelling  waves  in  the  Euler-Bernoul 1 i  beam  varies  inversely  with 
the  square  of  the-  wave  length,  "square  root"  damping  can  be  inter¬ 
preted,  regardless  of  the  boundary  conditions  in  force,  as  a  state¬ 
ment  that  the  damping  exponent  is  proportional  to  the  distance 
covered  per  unit  time  by  a  travelling  wave  of  the  particular  wave 
length  in  question.  However,  this  explanation  seems  rather  suspect 
to  the  writer  because  the  expression  of  standing  waves  in  terms  of 
travelling  wvee  is  a  bit  of  a  mathematical  artifice  here.  The  expla¬ 
nation  would  become  quite  appealing,  however,  if  it  could  be  demon¬ 
strated  that  in  long  thin  beams  the  attenuation  of  travelling  waves 
is  specifically  a  function  of  distance  travelled  rather  than  other 
factors .  The  author  is,  at  the  present  time,  unaware  of  any  research 
which  may  have  been  carried  out  in  this  direction. 


. v  guiding  philosophy  in  our  modelling  work  is  that  a  proposed 
mathematical  model  for  a  physical  system  must: 

(i)  Be  well-posed  from  the  mathematical  point  of  view; 

(ii)  Replicate  the  physical  phenomena  observed; 

(ii.it  Involve  only  equations  all  of  whose  terms  can  be  assigned 
a  direct  physical  meaning  in  terms  of  observable  system 
properties;  or  characteristics. 


Where  the  boundary  conditions  on  the  beam  are  not  paired,  the  model 
(2.3),  at  least  as  far  as  we  understand  it  now,  does  not  meet  the 
third  criterion  stated  and  must,  consequently,  not  be  accepted  as  a 
viable  model  at  the  present  time.  The  apparent  necessity  of  discar¬ 
ding  this  model  under  the  circumstances  just  described  represents  a 
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real  loss  because  the  system  (2.3)  has  very  attractive  mathematical 
properties,  as  outlined  in  [  E  ]. 

Since  the  middle  term  in  (2.7)  has  a  ready  physical  explana¬ 
tion  as  the  rate  of  change  of  bending,  it  is  tempting  to  try  to  use 
this  equation  whatever  its  disadvantages.  If  no  modifications  of 
the  boundary  conditions  are  made,  as  compared  with  the  Euler  -  Ber¬ 
noulli  conditions,  it  is  not  true  in  general  that  the  energy  is  mon¬ 
otone  decreasing.  Consider,  for  example,  the  cantilever  case  and 
the  initial  state 


w  ( x  ,  0  )  - 


0 ,  2= ( X , o ; 

at 


=  x2  -  -x3  +  -X4 


\  -c-  S  ( w  { 

•,t) , 

’,t)) 

1  =  2r  ( 

S  dt 

at 

r+ 

II 

O 

c.__ 

(Note  that  -0'-  is  only  required  to  lie  in  L  [0,L]  for  finite  energy 

at 

2  2 

solutions,  i.e.,  solutions  in  the  state  space  H  [0,L]  x  L  [0,L] 
appropriate  to  the  corresponding  semigroup,  and  thus  need  not  satis¬ 
fy  the  boundary  conditions  imposed  on  w  at  the  free  endpoint  x  =  L. 


9  w 

a  tax' 


so  that  the  energy  is  actually  increasing  during  some  interval  afte: 
t  =  0,  rather  than  decreasing. 

Undaunted  by  this,  it  is  tempting  next  to  see  if  some  modifi¬ 
cation  of  the  Euler  -  Bernoulli  boundary  conditions  will  result  in 
solutions  of  (2.7)  having  monotone  decreasing  energy  at  all  times. 
One  can  realize  this  by  modifying  the  zero  lateral  force  condition 
applying  at  an  end  point  to 


^2-w  -  El  -  *  =  0  . 

9 tax  9xJ 


(2.8) 


and  one  then  finds,  for  smooth  solutions,  that  the  rate  of  change 
of  energy  is  equal  to 


Unfortunately  our  satisfaction  is  short  -  lived,  for  we  realize  that 
the  beam  with  free  endpoints  has  purely  inertial  solutions 


w(x,t)  =  Wq  +  W^t  +  W2X  +  Wgtx 

o 

for  which  (2.9)  has  the  value  -2yLw3  <  0  for  w3  /  0  even  though  no 

energy  dissipation  should  take  place  in  such  motions;  we  have  over 
done  our  attempt  to  introduce  damping  here.  Nevertheless,  the  exer¬ 
cise  is  not  quite  a  total  loss.  If  we  work  a  little  harder  and 
change  (2.8)  to 


-El 

-3 
o  w 

o 

+  2r 

a3w 

-  li  f 

3" ( L , t ) 

_9w/ 

0,  t) 

]  =  0 

/ 

(2.10 

ax3 

at  ax 

T  l 

XJ 

St 

at 

J 

we  then 

f  ind , 

for 

smoo 

th  solutions 

,  that 

—  S(w( 

’  - 1 )  , 

9W(  . 

-t)  ) 

=  ~2  >■ 

[L(  a'W 

(x,  t )  j 

2 

dx  + 

2 

—  (  L  ,  t ) 

_£hj 

;o,t)l 

dt 

at 

t 

l0L  axax 

Ll 

■at 

at 

j 

and  the  energy  is  seen  to  be  non-increasing  for  all  solutions  while 
remaining  constant  for  the  purely  inertial  motions.  Condition  (2.10) 
seems  somewhat  unnatural  since  it  is  non-local.  In  fact,  we  do  not 
propose  (2.  1C)  as  a  completely  serious  solution  to  our  problem,  but 
this  small  success  leads  in  a  natural  way  to  the  more  believable 
model  which  we  proceed  to  discuss  in  the  next  section. 


> 
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3  .  Form  of  _t he  Mathematical  Model . 

As  we  have  noted  in  the  previous  section,  the  commonly  used 
Euler  -  Bernoulli  beam  equation,  which  will  serve  as  the  starting 
point  for  our  discussion  here,  is  based  on  conservation  in  time  of 
the  energy  expression  (1.1),  i.e., 


£  1  w ,  - - 

at 


[  M  *  )2  *  -  (  &  )2  ] 


wherein  w(x,t)  is  abbreviated  to  w,  p(x)  and  EI(x)  to  p  and  El,  res¬ 
pectively.  We  will  assume  here  that  p  and  El  are  uniformly  positive 
on  0  <  x  i  L ,  that  p  is  continuous,  at  least,  and  that  El  is  at 
least  piecewise  twice  continuously  differentiable  on  [0,L]. 

Assuming  for  the  moment  that  the  function  w(x,t)  describing 
the  evolution  cf  the  beam  displacement  is  smooth,  an  easy  calcula¬ 
tion  shows  that 
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The  presence  of  the  angular  velocity  expression  in  the  under- 

atax 

lined  term  indicates  that  its  coefficient,  -  —  f  El  ^  ]  ,  should 

ax  1  ax2  J 

be  interpreted  as  a  restoring  torque ,  arising  due  to  spatially  vari¬ 
able  bending  of  the  beam.  Realizing  that  this  coefficient  represent; 
a  torque  aids  us  in  interpretation  cf  th  damping  term  which  we  new 


*-  *.  %  "n.  •  ", 


introduce  into  the  system  via  the  definition 


rL 

x ,  t )  =  2  h  ( x ,  i  ) 

f  (X,t) 

-  ^  U,t) 

|  d*. 

(3.3 

Jo 

L  3t3x 

3t9x 

We  think  of  as  a  torque  acting  on  the  beam  at  the  point  x  due 

to  the  differential  rotation,  as  compared  with  the  rotation  at  x,  of 
tiie  beam  at  points  %  "near"  x.  In  many  cases  the  support  of  the 
" i nter rac t ion  kernel"  h(x,£)  would  be  restricted  to  a  thin  strip  in 

R" ,  centered  on  the  line  x  =  f.  Application  of  Newton's  second  law 
dictates  the  symmetry  condition 

h U,x)  =  h(x,$)  .  (3.4) 

In  constant  coefficient  applications  it  is  convenient  to  replace 
h(::,$)  by  yh(x-£},  where  y  >  0  is  used  to  parametrize  the  strength 
of  the  damping  effect  and  h(q)  satisfies  the  normalization  condition 


]  h( n)  dn  =  i  (3.5) 

0  _  c 

C 

and  the  even-ness  condition  h(rj  =  h(-q). 

The  source  of  the-  damping  torque  in  differential  rotation  is 
best  illustrated  for  the  case  of  beams  composed  of  composite  mater¬ 
ials,  such  as  fiberglass,  boron  and  graphite  composites  and  wood. 

We  may  imagine  that  long  fibers,  whose  modulus  of  elasticity,  per 
unit  cross  -  sectional  area,  is  greater  than  that  of  the  beam  as  a 
whole,  pass  through  the  beam,  held  in  place  by  a  matrix  material  of 
some  sort.  As  the  beam  undergoes  deformation  of  various  sorts,  beam 
elements  at  x  and  £  may  rotate  at  different  rates,  reflected  by  dif- 

ferent  values  of  {>:,  t),  - - (  £  ,  t )  .  if  we  think  of  the  fibers 

B 1 3x  B t  Bx 

themselves  as  having  nearly  constant  length,  differential  rotation 
must  result  in  movement  of  the  fibers  relative  to  the  matrix,  with 
accompanying  friction  against,  or  deformation  (largely  inelastic)  of 
the  matrix  material.  The  result  is  a  torque  of  the  type  which  we 
have  just  described,  since  the  net  motion  of  the  fibers  relative  to 
the  matrix,  within  the  individual  beam  elements,  will  be  different 
or:  one  side  of  the  elastic  axis  from  what  it  is  on  the  other  side- 


when  differential  rotation  is  taking  place  (see  Figure  3.1). 


Fig .  3 . 1 

We  add  the  terir. 

Tv  (  X  ,  t )  — — —  (  X  <  t )  dx 

n  dtdx 

to  both  sides  of  (3.2) ,  so  that 
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(3.6) 


Equating  the  separate  parts  of  (3.6)  to  zero  (a  procedure  which  can 
be  justified  by  the  principal  of  virtual  work)  yields  the  integro- 
partial  differential  equation 


a2w 

p  at2 


+  a_  [  a  r  EIa_w  -i  _  T 

ax  L  ax  L  j 
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a2w 

'"at2 
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Latax  atax 
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(3.7) 


and  the  requirement  that,  at  x  =  0  and  x  =  L, 


?  O 

EIa  w  a..w_ 

ax2  atax 


_  [  a_  f 


EI 


a2w 


ax 


)  - 
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Various  beam  configurations  now  lead  to  different  sets  of  boundary 
conditions.  For  example,  in  the  case  where  the  beam  is  clamped  at 
x  =  0  and  free  at  x  =  L  (i.e.,  the  cantilever  case)  we  obtain 


w(  0  ,  t  ) 


2 

0,  *w  (0,t)  =  0,  (L,t)  =  0, 

ax  ax^ 


(3.9) 


(l, t)  ' 


at-; 


V-  ( I. .  t )  t ,  t 


atax 


)Jd$ 


=  0 .  (3.10) 


The  reader  will  be  able  to  generate  boundary  conditions  correspond¬ 
ing  to  other  configurations  with  equal  ease. 

Equation  (3.6)  now  becomes 


(since  the  roles  of  x  and  £  are  symmetric  and  h(x,t)  =  h(£,x)) 


-  «  ■>  »  k  -w* 
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h  ( x ,  t )  F  5_w_  ( x  ( t )  -  1  — — — —  ( x ,  t )  d£dx 

L  9t9x  9 1 9x  -1  9t9x 


2  2  2 

h(x,af  -  -§-^U,t)l  d£  dx  <  0  .  (3 .11) 

1  9t9x  atdv  -1 


9t9x 


Moreover,  assuming  that  h(x,£)  is  not  identically  zero,  the  inequal- 

^2. , 

ity  is  strict  except  in  precisely  the  cases  where  — -  is  constant, 

9t9x 

i.e.,  the  inertial  motions.  Thus  energy  is  strictly  decreasing 
whenever  the  bending  rate  is  not  identically  zero  and  is  conserved 

Q 

when  the  bending  rate,  - - *,  vanishes  identically. 

9t9x^ 

In  general  the  form  of  h(x,£)  will  depend  on  the  elasticity  of 
the  fibers  of  the  material,  as  discussed  earlier,  the  distribution 
of  their  lengths,  the  nature  of  the  interface  between  them  and  the 
matrix  material,  etc.  The  more  elastic  the  fibers  and  the  shorter 
their  average  length,  the  more  the  "mass"  of  h(x,£)  will  be  concen¬ 
trated  near  the  line  $  =  x. 

For  perfectly  inelastic  fibers,  all  having  length  L,  the  same 
length  as  the  beam  itself,  h(x,£)  is  a  constant;  call  its  value  y/L. 
Then 
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h(x,£)  (£.t)d£  =  x  —  (L.t) 

atax  l  1  at 


(0,t) 


h(x,£)  — — — ■ — ( x , t )  d£  -  y  — w- ( x , t )  , 

atax  atax 


so  that  (3.7)  becomes 


_  2>  •'«  + 

at2  atax2 


2w  1  _ 

32  j  ' 


and  the  fourth  boundary  condition  of  the  cantilever  case  in  (3.10), 
for  example,  becomes 


t  EI(^ 


a2w  I 


+  2r  ( L ,  t )  -  2*  (L,t)  =  o  , 

x  =  l  atax  l  at 


in  exact  agreement  with  the  model  obtained  "ad  hoc"  at  the  end  of  th 
previous  section. 


In  the  case  of  metallic  beams  the  "fiber"  explanation  used 
above  is  not  persuasively  valid.  It  may  be  possible  to  think  of 
elongated  crystal  structures  in  the  material  playing  much  the  same 
role  as  the  fibers  above  but  the  analogy  may  well  be  far-fetched. 

The  evidence  for  frequency-proportional  damping  is  not  quite  as  comp 
elling  at  this  writing  for  metallic  beams  as  it  is  for  composite  and 
wooden  beams,  as  we  will  see  in  our  discussion  of  experimental 
results  in  Section  9. 


4 .  Additional  Properties  of  the  Model . 

Application  to  metallic  beams,  if  valid  at  all,  would  necess¬ 
arily  involve  the  supposition  that  the  fiber/crystal  structures,  in 
terns  of  which  the  damping  action  has  been  explained,  are  rather 
short.  It  is  therefore  of  some  interest  to  investigate  the  limitin 
form  of  the  equations  as  the  support  of  h(x,£)  is  restricted  to 
small  neighborhoods  of  the  line  x  =  $. 

To  avoid  complications  with  the  boundary  conditions  here,  let 
us  consider  an  infinite  beam,  so  that  0  is  replaced  by  and  L  by 
.  Then  the  energy  loss  term  in  (3.11)  is 

2  2  2 

h(x,$)[  (x,t)  -  dx  .  (4.1) 

L  3t3x  9t3x 

Using  the  change  of  variables 

y  =  -±z  U  +  t),  n  =  -4r  (X  -  £)  , 

J  z  </  z 

the  integral  (4.1)  may  be  rewritten  as 
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Now,  assuming  w  sufficiently  smooth, 
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and  the  integral  is 


2  [  f  h  (-  1  _  (  y+  h )  <  ( y-  H )  j  "2  [  [  “—[-==(  Y±on)  .t]  do  j2dr1dy.  (4.2 

v  z  2  '  J  n9t3x  2 


Consider  a  family  of  functions 


h  ( x ,  £  )  =  hy(x,i)  ,  y  >  0  , 

each  having  the  properties  assumed  earlier  for  h(x,£)  and  such  tha 
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where'  K(y)  and  G(y)  are  positive  integrable  functions  of  y  and 

1 im  J ( > )  =  0 . 

r-c 
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Assuring  that  ’  -  is  such  that 

9t3::2 


0^  ~3 


[  i  ]  =b 


on  C  V-'-  v 


uniformly  for  all  y  and  r;  while 
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uniformly  for  all  y  in  (-•,■)  and  for  rt  in  [—>,>]  ,  where 


1  i  m  £  (  y  )  =  0  , 
y->  0 


the  limiting  value,  as  >■ 


0,  of  the  integral  (4.2)  is  seen  to  be 


",  res- 


Returning  to  equation  (3,1)  with  0,L  replaced  by  -<*>, 
pectively,  we  see  that,  if  the  term  (4.3)  is  added  to  both  sides  of 
that  equation,  there  results 


— -  S  ( w  (  '  ,  t  )  ,  —  (  ’  ,  t )  )  +  2|  H(x)  f9-3^--(x,t)]2dx  = 
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+  EI  1 
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a  tax2 
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=  (after  integrating  the  last  product  twice  by  parts  w.  resp.  to  x ) 


r  H  siw —  i  +  r  E:  r_w  ]  ]  dx  . 

By.2  1  at  ax2  J  ax2  L  3x2  J  -1 

(Here  we  have,  of  course,  assumed  that  w,  —  and  their  x-derivative 

9t 

tend  to  zero  at  an  appropriate  rate  as  |x|  ->  « .  )  Now  the  condition 
that  the  integrand  should  vanish  identically  corresponds  to  the  Kel 
vir.  -  Voigt  partial  differential  equation 


a2w  +  2  a2 

f  K  2!-w-  -  1 

+  32.. 

j  EI  } 

=  0  . 

(4.4) 

at2  ax2 

1  at  ax2  J 

ax2 

L  ax2  J 

Thus  our  model  is  consistent  with  the  Kelvin  -  Voigt  model,  under 
the  stated  conditions,  as  the  "mass"  of  h(x,$)  is  progressively  nor 
concentrated  near  x  =  .  For  fixed  y ,  equivalently  fixed  h(x,£), 
this  means  that  our  model  approximates  the  Kelvin  -  Voigt  model 
for  vibrations  of  large  wave  length,  i.e.,  low  frequency.  This  wil 
be  made  more  precise  in  the  following  material  on  the  spectral  prop 
erties  of  the  constant  coefficient  equation  on 

Further  work  of  a  more  precise  character  will  be  required  in 
order  to  obtain  a  rigorous  argument  showing  that  solutions  of  (3.7) 
approach  solutions  of  (4.4)  as  y  -*  0  and  to  account  for  the  same 
phenomena  in  finite  beams  with  their  attendant  boundary  conditions. 
We  have  some  preliminary  developments  in  this  direction,  but  the 
work  must  appear  elsewhere. 
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We  pass  now  to  a  discussion  of  the  spectral  properties  of  the 
constant  coefficient  equation  on  (-”,■*').  Taking  our  earlier  remarks 
into  account ,  the  constant  coefficient  version  of  our  equation,  with 
r,  =  >:-£,  if; 


-  2} 


h(n)f  (K,t)  -  { x-  r> ,  t )  ]  d  t\  +  EI^~  =  0.  (4.5) 

L  9 tax  3ta>:  J  3x* 


ata>: 


The  function  h  has  the  properties  described  in  (3.4),  (3.5)  and  p, 

El  are  positive  constants.  Equation  (4.5  )  applies  for  an  infinite 
beam  or,  in  much  the  same  form,  for  the  periodic  case  wherein  x  is 
identified  with  x  +  I,  for  some  finite  L  >  0.  Our  purpose  in  this 
paragraph  is  to  study  spectral  properties  of  this  equation  under  the 

assumption  that  h  €  I, :  ( -  ,  )  with  h(r,)  dr;  =  1  .  Accordingly,  our 

J  _  . 

equation  becomes 
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h  ( n ) 


- — (x-n,t)  dr,  +  El  -  = 

3..  2  3..  * 


-<*-*  =  o,  (4.6 

ax4 


the;  form  which  we  will  use  for  our  analysis.  The  periodic  case  may 
be  developed  along  quite  similar  lines. 

We  begin  by  looking  for  solutions  of  (4.6)  in  the  form 

w ( x ,  t )  =  eot'f-(x) 

leading  to  the  steady  state  equation 


pc2'*  ( y- 


to  I  I  t  *-  ' 

t  \  / 


+  c  [  h(  rj<f'  '  (x-rjdq  +  EI<f(lv)(x) 


=  0.  (4.7 


e  1  ( x  )  dx, 


<  T  <  •>  , 


(4.6) 


be  the  Fourier  transform  of  <f ,  convergent  in  an  appropriate  sense. 
Applying  the  Fourier  operator  to  the  equation  we  have 


[  pc2  +  2  }■  t  2  (  1  -  h  {  t  )  )  c  +  EIt4  ]  +(t)  =  0.  (4.9) 

The  coefficient  of  + ( r )  must  vanish  on  the  support  of  the  function 
(or  distribution)  T(r),  leading  to  the  following  relationship  be¬ 
tween  <5  and  t: 


c  =  c  (  T  ) 


g(T)  +  J  g(r)2  -  pEl 


(4.10) 


with 


g  (  t  )  =  r  ( i  -  h  ( r ) ) . 


(4.11) 


Here,  as.  implied  earlier,  h(m)  is  the  Fourier  transform  of  the  int- 
erraction  kernel  h(rt).  From  the  integrabi  lity  of  h  ( rq )  we  know  that 

h(-r)  is  continuous  and 

1 im  h ( r )  =  0 .  (4.12) 

I  rh" 

For  a  relatively  low  level  of  damping  we  may  assume  y  <  pEI .  Then 
from  (4.10),  (4.11),  (4.12),  it  is  clear  that  for  large  values  of  t 
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to  first  order  in  h(m)  as  |t|  ■»  Asymptotically,  as  |t|  - 

r.  lies  along  the  rays 

R+  _  =  {  z  |  z  =  r(-r  1  iv'/  pEI  -  V2  ),  r  >  0  } 


in  the  left  half  complex  plane. 


From  (3.5)  it  is  clear  that 
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so  that 
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the  origin, 
show  that 


the  spectral  curve  is  tangent  to  the  imaginary  axis  at 
In  fact,  since  h  is  even,  two  integrations  by  parts 
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Thus,  at  longer  wave  lengths,  the  spectrum  tends  to  that  of  the 
corresponding  Kelvin  -  Voigt  equation,  as  we  should  expect  from 
our  earlier  work  in  this  section  since  taking  long  wave  lengths 
for  fixed  h  may  be  seen  to  be  equivalent  to  allowing  the  support 
of  h  (  r\ )  to  contract  to  (0),  keeping  wave  length  fixed.  A  typical 
spectral  curve  for  the  beam  of  infinite  length  with  the  proposed 
damping  mechanism  is  shown  in  Figure  4.1. 


For  finite  beams  of  large  length,  L,  the  spectrum  will  consist  of 

discrete  points  closely  spaced  along  a  similar  curve  in  the  complex 
plane. 

From  (4.9  )  we  have 

h(T)  =  i  +  -SHi  +  . 

2yc[r)  2yrz 

Thus,  in  principle,  h(r),  and  hence  h(r|),  can  be  reconstructed  from 
the  knowledge  of  c(t).  This  remains  true  in  an  approximate  sense 
for  finite  beams  of  large  length  and  may  prove  useful  in  identifying 
actual  dissipation  mechanisms  in  the  laboratory. 

It  is  interesting  to  note  what  happens  when  the  basic,  conserv¬ 
ative,  system  (before  damping  is  introduced),  i.e.,  the  Euler  -  Ber¬ 
noulli  model  in  the  above  analysis,  is  replaced  by  the  Rayleigh  mod¬ 
el  (1.5).  The  equation  (4.9)  is  then  replaced  by 

[( P  +  IpT2)c2  +  2yr2 ( 1 -h ( t) ) c  +  EIt4]*(t)  =  0 


which  leads  to  (cf.  (4.10)) 
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5 .  State  Spare,  Boundary  Conditions, _ Semigroup  Formulation . 

From,  the  form  (1.1)  of  the  energy  expression  for  the  Euler- 
Bernoulli  beam  we  are  led  to  define  the  state  space  as  a  subspace, 

2  2 

or  in  some  cases  a  quotient  space,  of  H  [0,LJ  x  L  [0,L].  The  inner 
product  should  be  consistent  with  the  bilinear  form 


< ( W , v ) , ( w , v ) >  =  J  p(x)v(x)v(x)  +  EI(x)w' ' (x)w' ' (x) Idx 

>1  n  L  J 


which,  in  turn,  is  related  to  the  energy-motivated  seminorm 


C(w, v; 


[p(x)jv(x)|2  +  El  {>:)  ;w"  (x)  !  2ldx.  (5.1) 


To  do  this  properly  it  is  necessary  to  say  something  about  the  boun 
dary  conditions,  not  all  of  which  are  created  equal.  Let  B  be  the 
set  of  distribution  pairs,  in  Schwartz's  sense,  such  that  the  boun¬ 
dary  conditions  on  the  unforced  beam  are  expressible  in  the  form 

,w>  •+  <$2'v>  =  0  «  e  B  '  (5.2) 

Some  of  these  boundary  conditions  represent  kinematic  constraints 
while  ethers  are  really  dynamical  equations.  The  former  constitute 
a  subset  C  C  D  of  pairs  (p,0),  where  p  is  a  continuous  linear  func- 

tional  on  H  [0,L],  which  prescribe  certain  constraints  on  the  beam 
displacement  without  any  reference  to  applied  forces,  such  as,  e.g. 


w  (  0  ,  t )  =  0,  32(0, 


=  0 ,  etc 


Let  us  denote  the  subspace  of  H^[C,L]  consisting  of  w  such  that 

2 

<p,w>  -  C,  (p,0)  C,  by  Hq  [  0 ,  L  ]  .  The  remaining  boundary  condi¬ 
tions,  corresponding  to  distribution  pairs  in  B  -  C,  are  really 
dynamical  equations  which  have  to  do  with  certain  forces  applie 
to  the  system  or,  in  the  free  case,  state  that  certain  forces  a: 
zero.  Thus,  in  the  cantilever  case 


are  statements  that  the  applied  torque  and  lateral  force,  respec¬ 
tively,  at  x  =  L,  are  zero.  In  the  forced  case  the  right  hand  sides 
would  be  replaced  by  a  moment  ;i  and  a  lateral  force  9.  Although 
the  elements  of  B  -  C  in  this  cantilever  case  still  have  the  form 
(  p  ,  0  )  ,  it  is  clear  that  p  is  not  a  continuous  linear  functional  on 

H2 [ 0 , L ] . 

2  2 

Let  Nq  be  the  subspace  of  H£[0,L]  x  L  [0,L]  on  which  the  semi¬ 
norm  (5.1)  vanishes.  In  many  cases  Nc  =  {0}  but  not,  e.g.,  in  the 

case  of  the  beam  with  two  free  endpoints,  or  in  the  case  where  one 
endpr int  is  free  and  the  other  hinged.  Extended  in  the  obvious  way 


f  [  H2[0,L]  x  L2[0,L]  ] //N( 


the  semi -norm  (5.1)  becomes  the  norm  for  our  state  space,  H.  We 
will  continue  to  refer  to  state  vectors  as  displacement  and  velocity 
pairs  ( w , v ) ;  the  equivalence  class  will  be  understood  in  those  cases 
where  Nc  ?  { 0 ) . 


The  state  space  now  defined,  we  proceed  to  pose  our  dynamical 
equations  in  that  content.  With  z  =  (w,v),  consider  the  unbounded 
operator  on  P  defined  by 


*.z  -  .  0  1  1 1  w  1 

l-fi  -G  J  l  V 


where  the  operational  definitions  c.f  A  and  B  are 


(  Aw  =  (  E  i  (  x  )  w  ’’(>:)  ) 


(Gt 


2'^  ;  h(x,U[v’(>:)-v’(t)]dC 


h'.e  domain  of  ft  consists  of  those  pairs  (w,v)  e  H4[0,L]  x  H2[0,Lj 


s:  h.  that  {5.7)  is  true  for  all  { Pj  ,  B .  Passing  to  equi¬ 

valence  classes  where  necessary,  this  is  a  dense  subspace  D  (  ft )  C  H 


and  A,  so  defined,  is  a  closed  operator.  Our  model  equation  (3.7) 
together  with  the  boundary  conditions  arising  from  (3.8),  and  dis¬ 
cussed  earlier,  is  then  equivalent  to 


z 


Az  . 


(5.3) 


The  computation  (3.11)  can  then  be  used,  with  a  few  refinements,  to 
show  that  A  is  dissipative,  i.e.. 


< z , Az>  <  0,  z  e  E( A)  . 

From  this  it  follows  immediately  (see,  e.g.,  [  F  ],  [  G  ])  that  A 

generates  a  strongly  continuous  semigroup,  Z(t),  of  bounded  operators 
or.  h  . 


It  turns  out  that  the  semigroup  Z(t)  is,  in  fact,  a  holomorp'nic 
semigroup,  as  one  might  expect  from  the  location  of  the  spectrum  in 
the  constant  coefficient,  infinite  length  case  discussed  earlier.  A 
general  idea  of  the  proof  of  holomorphicity  proceeds  in  the  follow¬ 
ing  way.  First  of  ail  one  studies  the  operator 


A 


o 


'  0  I  ' 
t-A  -Oj 


where 


( GQv)  (x) 


H ( x ) v 1 (x) 


H  ( x ) 


rL 

h ( x , t )  d£  . 

^  0 


the  domain  of  A  being  defined  in  terms  of  a  suitably  adjusted  set 
of  boundary  functionals,  B  .  The  boundary  conditions  corresponding 
to  the  subset  C  C  B  remain  unchanged  as  a  subset  of  B  ,  so  the  state 
space  Y.  remains  unchanged.  The  problem  of  finding  and  estimating  the 
normalized  eigenfunction  pairs  (#•  ,  i'^)  and  associated  eigenvalues, 

A^  ,  for  A  is  a  complicated,  but  entirely  standard  [  K  ],  eigenvalue 

problem  for  a  system  of  ordinary  differential  equations  with  bound¬ 
ary  conditions  of  standard  typ  „■ .  It  can  be  seen,  and  the  details 
will  appear  elsewhere,  that  the  A^  asymptotically  lie  along  the  rays 


arg(z)  =  +  (n/2  +  e) 

in  the  left  half  complex  plane,  where  e,  satisfying  0  <  £  <  n/2, 
depends  upon  p(x)  ,  H(x)  and  EI(x)  only,  not  upon  the  boundary  con¬ 
ditions  specified  via  BQ .  The  eigenfunction  pairs  (#^,4^)  may  be 

shown  to  form  a  Riesz  basis  for  ff.  It  is  possible  that  some  finite 
number  of  the  have  positive  real  part  because  A0  is  not  dissipa¬ 
tive  in  general.  The  Riesz  basis  property  of  the  ( §  ^ )  together 
with  the  asymptotic  location  of  the  are  all  that  are  required  to 
show  that  A0  generates  a  holomorphic  semigroup  on  H  (cf.  [  F  ]). 

Now  let  and  /p  denote  the  normalized  eigenf unction 

pairs  and  associated  eigenvalues  for  A.  It  turns  out  that  one  can 
prove  that 

lim  (  |  A,.  -  Xk|  )  =  0, 
i>.ki-  k  k 


From  the  last  relation,  using  a  standard  Hilbert  space  result  (see 
[  H  ],  p.  208  ff.)  it  can  be  seer,  that  the  (<Pj.,4j,)  also  form  a  Riesz 

basis  for  Y  so  that  A  is  seen  also  to  generate  a  holomorphic  semi¬ 
group  on  h . 


6 .  Spline  [  Fin:te  Element  Approximation . 

According  to  the  principle  of  virtual  work  [  I  ],  given  w(t), 
the  trajectory  of  a  mechanical  system  subject  to  holonomic  constr¬ 
aints  K ( w ( t ) )  =  0,  with  forces  f^(t),  i  =  1,2,  ... ,r,  acting  on  the 

system,  the  work  done  by  these  forces  in  a  small  virtual  displace¬ 
ment  from  w(t),  at  any  instant  t,  vanishes.  More  precisely,  if  v 
is  a  vector  such  that 

?K  ( w  ( t )  )  v  =  0 


then 


2  fi(t)  v  =  o 
i  =  l 


(6.1) 


It  should  be  noted  that  the  sum  (6.1)  must  include  all  forces  act¬ 
ing  on  the  system:  inertial,  restoring,  damping,  exogenous,  etc.  If 
a  damping  force,  linearly  dependent  on  velocity, 

g ( t )  =  -  Gw ( t ) , 

where  G  is  a  non-negative  self-adjoint  operator  (matrix  in  the  fin¬ 
ite  dimensional  case)  is  introduced  into  an  otherwise  conservative 

system  whose  kinetic  energy  is  ^w(t)*Mw(t)  and  potential  energy  is 

Aw(t)‘Vw(t),  M  and  V  also  non-negative  self-adjoint,  M  positive 
definite,  (6.1)  may  be  expressed  as 


j^-Mw(t)  -  Gw(t)  -  Vw(t)j  v  = 


where  -Mw  is  the  inertial  force,  -Gw  if-  the  damping  force,  and  -Vw 

is  the  restoring  force  for  the  system.  Replacing  v  by  w  (which  mus 
be  an  admissible  value  for  v)  we  obtain 


\  -v 


%  S.  \ 
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—  (w(t),w(t))  =  -w(t)*Gw(t)  (6.3) 

at 


where  C  is  the  energy 

C ( w ( t ) , w ( t ) )  =  j fw ( t ) *Mw ( t )  +  w(t)*Vw(t)J  .  (6.4) 

Formally,  the  derivation  of  the  equation  of  motion, 

Mw ( t )  +  Gw ( t )  +  Vw ( t )  =  0 ,  (6.5) 

is  carried  cut  via  (6.2);  in  practice  it  is  usually  inferred  from 

(6.3)  after  a  quadratic  form  -w(t)*Gw(t)  expressing  the  energy  loss 
through  damping  is  hypothesized.  In  spatially  distributed  systems  the 
forms  in  which  Q,  V  appear  in  (6.2)  may  be  different  from  the  forms 
in  which  they  appear  in  (6.3),  (6.5)  due  to  transformation  via  integ¬ 

ration  by  parts  (use  of  the  divergence  theorem),  appearance  of  bound¬ 
ary  terms,  etc..  For  example,  (6.3)  might  be  expressed  as  a  quadra- 

tic  form  in  - - ,  but  application  of  (6.2)  requires  that  v  should  be 

5x3t 

an  admissible  velocity  -  not  a  partial  derivative  of  such  a  velocity, 
necessitating  a  transformation  by  integration  of  parts. 

To  implement  a  spline,  or  general  finite  element,  approximation 
procedure  for  the  damped  beam  equation  we  begin  with  a  finite  dimen¬ 
sional  vector  §  of  coordinate  functions 


f  = 


We  may,  optionally,  assume  that 

<)3,'fi>  =  0,  i  =  1,2 . n,  (j3,0)  e  C, 

so  that  the  9^  satisfy  the  kinematic  constraints  from  the  outset,  or 

these  constraints  may  be  imposed  in  a  later  step  to  be  described. 

For  beam  applications  the  basis  functions  9^,  i  =  l,2,...,n,  must  be 


(6.6) 


such  that,  with 


w(x,t)  =  v  wi(t)'Pi(x)  =  §*w(t) 
i  =  l 


(6.7) 


(£,t)  =  2  wi(t)<Pi(U  =  #u: 

2  =  1 


(6.8: 


atax 


(t.t)  =  2  wi(t)9i  (£)  =  ft '  (£)*w(t) 


(6.9) 


are  well-defined  quantities  at  points  £  where  external  lateral  forces 
and  torques  are  to  be  applied. 


Use  of  the  representation  (6.7),  with  concomitant  expression 
of  etc.,  as  in  (6.8),  (6.9)  implies  movement  of  the  beam  is  sub- 

at 

ject  to  certain  holonomic  constraints;  infinitely  many  of  them,  in 
principle,  the  net  effect  of  which  is  to  allow  the  system  to  move  in 
the  finite  dimensional  subspace  spanned  by  the  coordinate  functions 
,  i  =  l,2,...,n.  With  such  representation  of  w(x,t),  the  work  - 

energy  relation  with  external  lateral  force  c  and  torque  t  applied  ai 
the  point  x  =  £  may  be  seen  to  be 


— -fw(  t )  *Mv:(  t ) +w( t ) *Vw( t )1 -w( t ) *Gw( t )  =  cf(£)*w(t)+T§,(£)*w(t).(6.10) 
dt  ^  J 


The  functions  9^  are  often  chosen  to  be  independent  basis  functions 

but  this  is  not  always  convenient  at  the  outset.  If  the  9^  are  not 

independent  there  will  be  a  set  of  r  independent  linear  equations 
satisfied  by  the  accompanying  coefficients,  which  we  may  represent  in 


the  form 


P  w  ( t )  =  0  , 


6.11) 


where  P  is  an  r  x  n  matrix,  r  %.  n  -  m  <  n.  Let  the  m  x  n  matrix  Q 


be  such  that 


is  nonsingular  and  define  Q  w  =  z.  We  then  have 


( T , U ) ,  we  find  that  w(t)  =  U  z(t)  and,  with 


Writing  f~l  = 

K  f  U*MU,  V  =  U*VU,  G  =  U*  GU  , 

we  have,  in  place  of  (6.10), 

f  z  ( t ) *Mz ( t )  +  z ( t ) * Vz ( t )  ]  -  z(t)*Gz(t) 
dt  L  J 

=  c(t)  §U)*Uz(t)  +  r(t)  *'(t)*Uz(t)  ,  (6.12) 

leading  to  the  equations 

K  z(t)  +  G  z(t)  +  V  z(t)  =  c(t)  U*«(U  +  t  ( t )  U*f'U)  (6.13) 

as  the  equations  of  motion  for  the  system. 

The  kinematic  boundary  conditions,  corresponding  to  C  as  descri¬ 
bed  above,  may  be  imposed  via  (6.11)  if  desired.  The  other  boundary 
conditions,  corresponding  to  B  -  C,  are  never  imposed  explicitly;  they 
come  about  as  particular  cases  of  the  dynamical  equations  in  (6.13). 

If  the  coordinates  to  be  retained  are  a  subset  of  the  components 
of  w,  we  can  write  w*  =  (w*,  w*  )  ,  where  w  comprises  the  components  to 

be  discarded  and  w  the  components  to  be  retained.  Then  P  is  partit¬ 
ioned  accordingly:  P  =  (P,  P),  and  we  can  assume  that  the  r  x  r  matrix 
P  is  nonsingular.  Then,  with  m  -  n  -  r , 

w  =  -(P)_1P  w  ,  w  =  [-^PT  w  ~  U  w 

and  we  obtain  (6.13)  again,  replacing  z(t)  by  w(t).  Computationally 
this  is  easier  because  P  is  r  x  r  while  [q]  is  n  x  n. 

In  actual  computational  work  it  is  preferable  to  define  the 
energy  and  dissipation  forms  in  the  continuous  context  and  then 
approximate  them  by  quadratic  forms  on  finite  dimensional  spaces. 


y 


Application  of  (6.11),  (6.12)  to  the  approximate  energy  and  dissi¬ 

pation  forms  then  leads  to  the  finite  dimensional  approximating 

equations  of  motion  in  the  form  (6.13)  with  symmetric  matrices  M , 

G,  and  V.  This  procedure  is  much  easier  and  more  reliable  than 
first  deriving  the  equations  of  motion  as  partial  differential  equa 
tions  ir.  the  continuous  case  and  then  proceeding  to  solve  those 
equations  numerically,  taking  into  account  all  of  the  boundary  cond 
it  ions.  In  fact,  by  the  recommended  procedure,  only  the  kinematic 
boundary  conditions  are  actually  enforced  explicitly.  The  example 
of  the  next  section  should  be  helpful  in  clarifying  these  matters. 


'  V  *  -  *' 
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7 •  Approximation  and  Simulation  for  the  Cantilever  Beam  System. 

Let  the  interval  k  =  [0,L]  be  divided  into  m  subintervals 
Ik  -  I"  ( k- 1 )  b  ,  kbl  =  [xk_, ,xfc]  =  [>:k-1  ,xk_1+h]  ,  h  =  L/m,  k  =  1,2, 

and  consider  general  quadratic  functions  on  the  intervals  I,: 


<3k(x)  =  ak  +  bk(>:'xk-l)  +  ck(x'xk-l'  * 

o 

The  functions  defined  on  Ik  by  1,  x-xk_j,  (x-xk_j)  '  anc^  extended 
by  0  to  the  rest  of  k  correspond  to  the  functions  ,  i  =  l,2,...,n, 

n  =  3m,  discussed  in  the  previous  section.  With 

a ( t ) *  =  ( a  J  t ) . ar ( t )  ) ,  b ( t ) *  =  ( b i ,  . . . , bm ( t ) ) , 


C  ( t ) *  =  ( c, { t  c  '  t  )  )  , 


(>:)  =  !  X- 


Xt  , x’Xt  .. 
-m  A  i 


-Xt  -x2’Xi 


• ' (x~xm-l)£ 'XI  I* 


and  w(tj"  =  (a  '  t  )  "' ,  b  ( t )  11  ,  c  ( t )  *  )  ,  we  can  define  a  function  w(x,t), 
piecewise  quadratic  in  x,  with  possible  discontinuities  at  the 
points  ::k ,  k  -  l,2,...,n-l,  by 

w  ( x ,  t )  =  ?  ( x )  *  w  ( t )  . 

We  proceed  to  define  the  energy  and  dissipation  forms  in  terms  of 
these  functions,  ignoring  the  discontinuities  for  the  present. 

Assuming  the  beam  to  have  mass  density  p,  the  kinetic  energy  is 


p 1  o 

-  t  \  ;  1)1  <2.-:  - 

u  o  '  - r 


m  „x,. 

P  r-  I  ^ 

-  “  I 
k=li-  x. 


[a-.  ( t )  b-_  it)  ( x-xk_ :  )  +ck  ( t )  (  x-.-j._j  )  2]  d: 


2  1  h*k  ( t )  "  + '  bj.  ( 1 )  2  +  3  ’  ck  ( t )  +h‘ak  ( t )  bk  ( t )  +  §h 3ak  ( t }  ck  ( t )  +ih  * bk  ( t )  ck  ( t )  j 

=  w(  t )  *Mw (  t  )  , 


from  which  we  see  that 


f  hl„ 

1h2l 

i  u  3  r 

m 

2 

m 

3 

m 

1hzIrn 

-  h3 1 

d-h4  I 

z  m 

3  m 

i  m 

-1  h3 1„ 

-1  h4l 

-  h  5 

l  s  m 

4 

m 

5 

m  j 

He:  e  Iv  denotes  the  m.  x  m  identity  matrix. 

Next  we  pass  tc>  the  potential  energy  which,  for  constant  bending 
modulus  El,  is 


f :  •  *  '  * 

j  ?:-w 

/ 

4  \ 

t  l-  J 

r  ns  r 

1  d::  =  iEI  I  1 

-  hx2 

J  * 

k=i  l’x, 


4ck(t)‘dx 


=  2 ElhCj.  ( t ) 


so  thru 


0  0 


V  =  0  0  0 

[  0  0  4EIhIr 


Passing  now  to  the  dissipation  form,  we  first  of  all  take  the 
interract  ion  kernel  h(x,£)  to  have  the  form  >h(x-£)  where,  assuming 
that  m  will  always  be  a  multiple  of  fixed  positive  integer  p , 

h  (  n )  -  1,  |  |  <  L  /  p  ,  h(n)  ~  0  otherwise. 

We  define  ,  ~  ( r. )  =  m/p  for  all  values  of  m  used.  For  the  case  m  - 
9,  v  =  p  ~  3 ,  the  support  of  h  (>:-£)  is  the  shaded  region  shown  in 
Figure  7.1. 


IV 


Fig.  7.1 


k  *  {  <X.S)  |  xk-l  s  x  s  xk  ,  xl-l  S  «  S  xi  } 


Then  in  Si,  we  have 


— —  (x,t)  -  U  ,  t }  =  bk(t)+2ck(t)  (x-xk._1)-bi3(t)-2c£  (t)  (  ) 

3t?x  3t<9x 

For  |k-.f|  <  i<  we  square  this  quantity  and  integrate  over  Sj.;  then 
sum  over  the  indicated  values  of  k,i  to  obtain  the  quadratic  form 


w(t)*G  w(t)  =  2 >  I  (a(t)*,b{t)*,c(t)*)  0  0  0 

k ,  £  n  v.2  -r  h3J 

0  h  Jk(  £  h  Jk  ( jj 

•  0  h5jk^  h*Klc,  i? 


f  a  ( t ) 


where 


[Jk,i]k  -  [Jk ,  ij]  J  -  1(  [Jk ,  jc] 5  =  (Jk,n)k  "  1  ' 


-*t vy vT  -.'vy -  •  ' 


.  ■TV- 


.  /  .*  .■  /  _* 


0  otherwise; 


=  -1  , 


K.*)k  *  v  K.JI  *  i-  K,(]?  *  K,J'-  -  -i  ■ 

^Kk  f  |  j  =  0  otherwise  . 

Additionally  we  must  treat  the  upper  and  lower  borders  of  the  region 

of  support  of  h(x-U ,  consisting  of  triangles  Tk,  T^,  k  =1,2,...,M -v 

as  shown  in  Figure  7.1.  After  integration  of  the  same  function  as 
shown  earlier  over  these  regions,  followed  by  summation  over  the  rel 
evant  indices,  we  obtain 


w  (  t )  *  0  ,  w  (  t  )  = 


2>  v  (a(t)*,b(t)*,c(t)*) 

k=  1 


•  0 

0 

0  1 

j' 

0 

b2  Jk ,  £ 

h3Kk,£ 

.  0 

h‘K,.  . 

rv  ,  .t 

h+Lk , e  J 

L 

a  ( t )  1 
b(t)  ; 

c(t)  J 


with  J,.  ^  arid  Kj.  ^  as  above  arid 


':k,  fjk  “  2  *  [Lk  ,  J}]  f  "  2  ’ 

LLk,  f]  j  = 


-  l~k,£j£  " 

=  0  otherwise 


r  ^  0 
-1-  tLk,fjk 


After  G  and  G  have  been  computed  as  indicated  we  set  G  =  Gn  +  G,  . 

Kent  there  is  the  matter  of  the  construction  of  the  reduction 
matrix  U.  From  the  representation  on  Ik 


w !  >: ,  t )  =  ak(t)  +  bk(  t)  (x-xk_1  )  -1-  ck(  t )  (x-xk_1  ) 


!  x  ,  t  )  =  b.  (  t )  a  2  c.  ( t )  ( x-x,.  _ 


we  see  that  w ( x , t )  is  piecewise  C‘  on  i  just  in  case,  for  k  =  1,2, 
. . . , m- 1  we  have 
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a  V.  .  ^  ( -t )  =  a,,  (t)  4  b.  ( t )  h  4  C.  (t)h2  , 


bk4-l(t)  =  kj.(t)  4  2ck(t)h  . 

For  the  cantilever  case  of  our  example  we  have,  in  addition, 

a . ( t )  =  0,  b i { t )  =  0  . 

Accordingly,  we  have  the  following  linear  relationships  holding  be 
tween  the  coefficients 


c  - : 


0  0  0, 
0  -h  0, 
0  -h  0. 


0  0  0  c 

0  0  -h2  0 

oo  o  -h2 


0  0  0 


0  0  0 


0  0  0 


Let  P  be  the  matrix  with 


0  0 


■1  1  0 ...  0  0  - 2h  0  0 ...  0  0 

0  -1  1  ...  0  0  0  -2h  0 ...  0  0 

0  00. ..-11  0  0  0.. . -2h  0 


Rl  =  -1,  j  =  i  +  1,  ri  =  0  otherwise 


Then,  with  Q  =  I  4  R  we  have 


I  4  R  hR 


I  4  R 


Q  hR 
0  Q 


Ly  inverted,  giving 


I  -  a. 
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(  *  j 


Q  1  -h(Q  1  -  Q  2; 


Since 


[  0 

Q  1 

f h-2R  ] 

_  '  h2  (  Q-Ijp  )  • 

l  2nR  J 

-  2h(Q-Im)  J 

we  obtain,  as  the  reduction  matrix 

h2I 


U  = 


-P  "P 


m 


3h2°  1 


; Q  1  +  2h  Q 
2  h  (  Q  ~  1  - 

I. 


I_) 

m 


m 


further  easy  computation  shows  that 


1 


r  i 

0 

0  . 

•  0  ] 

■  1 

0 

0  . 

.  0  ' 

j- 

i 

0  . 

.  0 

2 

1 

0  . 

.  0 

-1 

: 

l 

-  • 

n-2  - 

'  g 

3 

2 

1 

.  0 

: 

- 

l 

1  . 

.  1 

m 

m-  1 

n-2  . 

•  1  . 

We  see  therefore  that  the  reduction  matrix  U  can  be  computed  expli 
itly  without  having  to  resort  tc  numerical  matrix  inversion  in  thi 

example.  Once  this  has  been  done  the  matrices  M,  G  and  V  of  (£.12 
can  be  calculated  readily  and,  using  a  Cholesky  decomposition  to 

write  K  in  the  form  Mf*,  with  £  lower  triangular,  the  equation  (6. 

can  be-  pre-  and  post -mul  tipled  by  and  (£*)-1.  respectively,  to 

yield  an  equa.tion  for  which  the  corresponding  first  order  system 


matrix  is 


a  = 


l  -I  :V(T*)  1 


V'  1  ~  !  V  +  \  i 

^  V  *-  / 


We  have  calculated 


mat  r  1:: ,  and  subsequer.  1 1  y 


uec ,  for  a  cantilever  beam  with  p  =  1,  y  =  .Cl  ,  El  = 

L  =  1C.  The  interraction  kernel  h(x~t)  was  taken  to  b 
supper 


its  ex 
10  0  a  r.  c 


restricted  to  (x-t;  =  L/2.  In  Figures  7.2 (a) 


wit  k 


ar.c 


we  show  the  eigenvalue  patterns  an  the  upper  left  half  plane  for  it. 
=14,  16  and  24  (not  all  of  the  eigenvalues  are  shown  in  the  last 
case).  Figure  7.2(b)  magnifies  the  lower  right  hand  corner  of  Figure 
7.2(a)  in  order  to  show  the  pattern  of  the  eigenvalues  corresponding 
to  low  frequencies  more  clearly.  Variational  considerations  account 
for  the  decrease  in  the  natural  frequencies  as  M  gets  larger  but  the 
increasing  negative  real  parts  of  eigenvalues,  or  exponential  damp¬ 
ing  rates,  came  as  somewhat  of  a  surprize.  This  property,  if  it 
can  be  substantiated  in  general,  would  be  quite  desirable  for  appli¬ 
cations  because  it  would  mean  that  damping  rates  are  systematically 
underestimated  in  approximations  of  this  sort,  providing  a  safety 
factor  in  control  applications.  From  Figure  7.2  it  is  seen  that  the 
quadratic  dependence  of  the  damping  rate  on  frequency  at  the  low  end 
of  the  spectrum  has  only  the  most  limited  range. 


8.  Some  Theory _ Re lated  to  Laboratory  Experimentation. 

The  foregoing  results  from  mathematical  simulation  can  be,  anc 
have  been,  compared  with  related  data  taken  from  elastic  beam  exper¬ 
iments  ir.  a  laboratory  setting.  This,  of  course,  is  nothing  new  by 
itself  but  perhaps  we  can  claim  that  we  are  in  a  position  to  look  at 
such  data  in  a  different  light  than  most  of  our  predecessors.  The 
author  is  fortunate  to  have  access  to  extensive  electronic  equipment 
for  collecting  and  analyzing  data  in  the  UW  MIPAC  (Modelling,  Infor¬ 
mation  Processing  and  Control)  Facility  and  great  advantage  has  been 
taken  of  this  experimental  capability  during  the  last  two  years. 
Before  describing  the  results  of  this  experimentation  program  and 
the  relationship  of  these  results  to  the  mathematical  simulation 
discusser,  above,  we  digress  briefly  to  describe  the  experimental 
context  and  some  of  the  inherent  problems  attendant  upon  the  obser¬ 
vation  of  vibrating  beams  in  the  laboratory. 

The  laboratory  beam  configurations  easiest  to  realize,  at  least 
approximately,  are  those  involving  either  clamped  ends  or  free  ends. 
The  results  which  we  cite  in  Section  9  were  obtained  in  the  clamped 
f re***  (cantilever)  setting  and  in  a  "pseudo"  free  /  free  config¬ 
uration.  The  clamped  /  clamped  beam  is  easy  enough  to  work  with, 
and  we  have  conducted  some  experiments  in  this  context,  but  we  do 
not  report  on  this  case  here,  mainly  because  we  have  only  one  clamp 
really  adequate  to  the  task.  That  the  task  of  realizing  a  good  ap¬ 
proximation  to  a  clamped  beam  endpoint  is  not  at  all  trivial  will  be 
commented  on  at  some  length  later  in  this  section.  The  free  /  free 
case  is  qualified  by  the  word  "pseudo”  because,  in  fact,  the  beam  is 
actually  suspended  vertically  by  a  very  fine  nylon  thread.  The  two 
experimental  conf icurat ions  are  illustrated  in  Ficure  8.1(a), (b;. 


Let  us  begin  by  discussing  this  pseudo  free  /  free  case;  from 
now  on  we  simply  refer  to  it  as  the  pseudo  free  configuration.  To 
obtain  a  reasonable  facsimile  of  the  free  beam,  whose  actual  reali¬ 
zation  requires  a  gravity-free  environment  achievable  only  in  space 
or  free  fall  conditions  with  accompanying  difficulties  as  regards 
excitation  and  measurement,  we  suspend  an  elastic  beam  from  a  suit¬ 
able  support,  preferably  a  high  one,  by  a  long,  very  fine,  nylon 
filament;  for  our  purposes  3  lb.  test  nylon  fishing  line  served  ad¬ 
mirably.  The  accelerometer  is  attached  to  the  beam  at  the  lower  - 
truly  free  -  end.  Just  attaching  the  accelerometer  introduces  a 
variety  of  problems  in  itself  but  we  avoid  that  difficulty  in  the 
present  article. 

We  proceed  now  to  discuss  the  coupled  beam  /  filament  system 
in  some  detail.  We  will  do  this  as  if  the  beam  were  supported  by 
a  single  filament  although  a  double  filament  arrangement  such  as  is 
shown  in  Fig.  8.1(b)  is  used  in  practice. 

We  are  interested  in  energy  losses  from  the  beam  through  the 
supporting  filament,  so  we  will  model  the  beam  itself  via  the  energy 


conservative  Euler-Bernoul 1 i  equation.  Assuming  the  beam  to  have 
length  L,  uniform  linear  mass  density  p  and  bending  modulus  El,  all 
constant ,  and  assuming  the  cord  supporting  the  beam  to  have  length  k 
and  linear  mass  density  r,  we  have  as  the  approximate  energy  expres¬ 
sion  for  the  s y s t em 
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with  the  constraint  w(L,t)  =  v(0,t).  Here  w  is  the  beam  deflection 
as  a  function  of  x  and  v  is  the  string  deflection  as  a  function  of 
s.  From  this  energy  form  one  readily  derives  the  equations 
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=  0,  0  <  x  g  L  , 


=  0,  0  <  s  <  k  , 


8.1) 


(8.2 


and  the  boundary  condi tons 
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At  s  = 


we  assume  no  reflection: 
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It  is  not  supposed  that  the  last  condition  is  completely  realistic, 
It  simply  replaces  more  complicated  assumptions  ensuring  that  very 
little,  if  any,  energy  entering  the  the  supporting  filament,  or 
string,  is  reflected  back  to  the  beam.  We  are  assuming  that  waves 


in  the  string  move  only  in  the  direction  from  s  =  0  to  s  =  & .  Then, 
in  fact,  equation  (8.5)  is  valid  with  $  replaced  by  s ,  0  £  s  £  £  ; 
in  particular  it  is  true  at  s  =  0.  Using  this  in  (8.4)  along  with 
the>  constraint  w(L,t)  =  v(0,t)  we  arrive  at 

El  ( I. ,  t )  +  j"  2P-  -  (rpL)l/2]  ( L ,  t )  =0  .  (8.6) 
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Since  we  assure  £  large  and  r  small, 
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is  a  small  positive  number  and  the  boundary  condition  (8.6)  becomes 
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If  we  define 


a  =  (  P  /  E I )  1  7  * 

and  assume  solutions  of  the  form 
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with  w  complex,  the  resulting  eigenfunction  equation  is 

w(lv)  -  (au)4  w  =  0  ,  (8.8) 

Applying  the'  boundary  conditions  (8.3)  and  (8.7)  to  the  general 
solution  of  (8.8),  and  setting 

X  =  a,-  ,  c  =  T/'a2  , 


we  obtain  the  determinantal  condition 
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which  simplifies  to 

2  A  -  2 A  cos  AL  cosh  AL  +  2d  cos  AL  sinh  Ah  -  2ci  sin  AL  cosh  Ah  =  C 
If  we  let 
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and  substitute  into  this  equation,  we  find  that,  to  first  order  in  8 
we  have 
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and  from  this  it  is  easy  to  see  that,  to  first  order  in  1/k, 
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Then,  to  zero  order  in  1/k, 
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Asymptotically,  then,  the  damping  of  the  k-th  mode  of  the  beam 
due  tc  the  presence  of  the  supporting  filament  is  uniform  and  the 
damping  exponent 


2c  =  2 Ely  =  2EI  f(rpL)i/a  _  S£L 
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tends  to  zero  as  r  a  o  and  £  a  <* ,  i.e.,  as  the  supporting  cord  be¬ 
comes  very  light  and  very  long.  We  judge,  therefore,  that  this  type 
of  support  mechanism  will  not  lead  to  damping  rates  comparable  at 
high  frequencies  to  those  arising  from  internal  sources  in  our  labor 
atory  experiments. 

A  somewhat  similar  analysis  shows  that  losses  due  to  excitation 
of  acoustic  waves  in  the  surrounding  atmosphere  by  the  high  f requeue 
vibrations  of  the  beam  will  also  be  asymptotical ly  uniform  and  thus 


not  likely  to  be  confused  with  internal  damping  effects. 


It  has  been  recognized  by  various  experimenters  for  some  time 
that  clamped  beam  experiments  are  subject  to  substantial  degrees  of 
interference  from  the  interraction  between  the  beam  and  the  nearly, 
but  not  quite,  rigid  clamping  device.  In  order  to  have  some  means 
for  evaluating  the  likely  degree  of  error  in  damping  experiments  un¬ 
der  these  circumstances,  we  ask  the  reader  to  consider  with  us  a 
clamping  arrangement  in  which  the  beam  to  be  tested  is,  as  shown  in 
Fig.  8.2,  the  relatively  thin  beam  lying  in  the  region  x  >  0,  ex¬ 
cept  for  a  portion  embedded  in  the  more  massive  beam  which  occupies 
the  region  x  <  0  and  which  serves  in  our  "thought  experiment"  (our 
apologies  to  the  memory  of  the  originator  of  that  useful  expression!) 
as  the-  clamping  device,  stiff  but  not  entirely  rigid  -  the  unintended 
result  of  most  of  our  laboratory  clamping  attempts. 


W®  denote  lateral  displacement  of  the  thin  beam  by  w(x,t),  x  ^  0, 
an-  lateral  displacements  of  the  thick  beam  by  W(x,t),  x  £  0.  We 
will  study  "monochromatic",  i.e.,  single  frequency,  solutions  of  the 
combined  system  under  the  assumption  that  "outgoing"  waves,  moving 
toward  -  ,  in  the  clamping  beam,  are  eventually  absorbed  and  not  ref¬ 
lected  back  to  x  =  0. 


Ado; ting  complex  notation,  we  represent  the  displacement  of  the 


thin  bear,  by 


w(x.t)  =  +  yei(u'2t+awX) 


(8.9) 


where  a’  =  (p/Et),  p  and  el  being  the  linear  mass  density  and  bending 
modulus,  respectively,  of  the  thin  beam.  The  first  term  represents 
waves  moving  outward  from  x  =  0  to  while  the  second  term  represents 
waves  moving  inward  toward  >:  =  0  from  Because  of  our  assumption 

that  all  waves  in  the  thick  beam  originate  at  x  =  0  and  are  not  ref¬ 
lected  back  subsequently,  we  represent  the  thick  beam  displacement  by 

W  ( x  ,  t )  =  be1  (  w^t+Ac;X)  (  (8.  10) 

where  A4  =  (E'EI),  E,  El  the  linear  mass  density  and  bending  modulus, 
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,  p  =  A  - 


EIA‘(A+a)  A  +  a 

t  is  convenient  to  renormalize  p  to  A  -  a,  y  to  A  +  a  so  that 
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If  the  two  beams  are  of  the  same  material,  so  that  E  =  e,  then,  since 
l  becomes  small  more  rapidly  than  p  as  the  thickness  of  the  thin  beam 
decreases,  we  shall  have  0  <  A  <  a,  -1  <  p  <  0  so  that  the  outgoing, 
i.e.,  reflected,  wave  in  the  thin  beam  has  uniformly  smaller  modulus 
than  the  incoming  wave,  independent  of  the  frequency  /  wavelength. 


Computing  the  boundary  values  of 
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with  the  formula  (8.11)  for  w(x,t),  we  arrive  at  dissipative  frequency 
dependent  boundary  conditions 
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It  may  be  computed  that  the  energy  associated  with  the  travelling 
wave  of  wave  length  2rr/c«,  call  it  C(o,t),  computed  over  a  suitable 
finite,  moving  domain  with  0  as  its  left  hand  boundary  in  order  to 
avoid  Z  being  infinite,  decays  according  to 

dH(u' =  -8Apw5 


while  the  energy  itself,  over  a  large  domain  as  described,  may  be 

seen  to  be  proportional  to  w4 .  We  conclude,  as  a  result,  that  the 
energy  decay  rate  exponent  must  be  directly  proportional  to  the 

square  root  of  the  frequency,  w2 ,  under  these  circumstances.  It  is 
clear  that  such  a  damping  relationship  may  seriously  interfere  with 
accurate  measurement  of  the  internal  exponential  damping  rate,  even 

though  we  expect  the  latter  to  be  proportional  to  u2  rather  than  w. 
While  actual  clamping  mechanisms  are  not  really  thick  beams  as  sup¬ 
posed  here,  this  analysis  is,  we  feel,  sufficiently  realistic  to  at 
least  partially  explain  the  very  severe  difficulties  encountered  in 
trying  to  measure  the  damping  rates  of  the  natural  modes  of  clamped 
beams  in  the  laboratory. 


In  a  laboratory  experiment  the  beam  as  excited  by  striking  it 
smartly  in  order  to  set  it  into  vibration.  A  record  such  as  the  one 
shown  in  Figure  9.1  is  obtained;  in  this  case  the  record  exhibited 
is  the  graph  of  an  input  voltage  to  UW  MIPAC’s  HP5451C  System  Analy¬ 
zer  proportional  to  the  free  endpoint  lateral  velocity  of  a  canti¬ 
lever  beam  clamped  at  the  other  end. 


Fig.  9.1 

The  record  initially  shows  a  substantial  high  frequency  component,  as 
evidenced  by  the  jagged  character  of  the  graph  on  the  left  hand  side. 
As  the  motion  continues  we  observe  that  the  curve  becomes  smoother, 
evidence  that  the  high  frequency  component  is  becoming  smaller  in  re¬ 
lation  to  the  overall  amplitude  of  the  motion.  In  fact,  the  amplitude 
of  the  fundamental  mode  is  not  very  different  at  the  end  of  the  time 
interval  from  what  it  was  at  the  beginning;  very  little  damping  is 
observed  at  this  low  frequency. 

To  obtain  a  better  appreciation  of  the  quantitative  factors  in¬ 
volved  in  the  evolution  of  the  beam’s  motion,  the  record  is  divided 
into  a  number  of  successive  segments  of  equal  length.  The  (Fast) 
Fourier  transform  of  each  segment  is  computed  and  multiplied  by  its 
conjugate  to  yield  the  power  spectrum  for  that  interval .  For  most 


purposes  it  is  desirable  to  form  the  logarithm  of  the  power  spectrum, 
which  we  henceforth  call  the  logarithmic  power  spectrum.  This  is 
what  is  displayed  in  Figure  9.2  (though  for  a  different  input  record 
than  the  one  shown  in  Figure  9.1).  Here  two  successive  logarithmic 
power  spectra  are  superimposed  to  yield  the  complete  diagram.  The 
peaks  evident  in  each  graph  correspond  to  the  frequencies  of  the 
natural  modes  of  vibration.  Since  it  is  the  logarithm  which  is  dis¬ 
played,  the  vertical  gaps  between  peaks  of  successive  spectra  are 
proportional  to  the  (negative)  real  part  of  the  complex  system  eigen¬ 
value  associated  with  that  mode  of  vibration  while  the  horizontal 
position  of  the  peak  corresponds  to  the  imaginary  part  of  the  eigen¬ 
value,  i.e.,  to  the  frequency  of  that  mode. 


Fig.  9.2 

Figures  9.2  through  9.4  display  logarithmic  power  spectra  from 
a  "pseudo-free"  boron-epoxy  composite  beam,  a  clamped  wooden  beam, 
and  a  "pseudo  -  free"  metallic  beam  ,  respectively.  In  these  fig¬ 
ures  we  have  added  short  horizontal  bars  to  more  clearly  indicate  the 
successive  peak  positions.  The  data  is  averaged  over  many  experiment 
using  the  same  beam,  recording  data  over  an  interval  of  the  same 
length,  but  exciting  the  beam  in  different  ways  in  order  to  adequate¬ 
ly  sample  all  modes.  The  approximately  linear  relationship  between 
damping  ratio  and  frequency  is  strikingly  evident  for  the  wooden  and 
composite  beams,  less  so,  but  still  fairly  convincing,  for  the  metal¬ 
lic  beam. 


1 


HZ  LIN 


Fir.  S.3 


Fig.  9.4 


It  might  be  natural  as  a  result  of  this  evidence  of  frequency 
proportional  damping  in  elastic  beams  to  suppose  that  all  comparable 
elastic  structures  exhibit  the  same  sort  of  relationship  between 
damping  and  frequency.  Indeed,  this  was  the  author's  initial  predis¬ 
position.  It  is  here  that  continued  experimentation  in  the  laboratory 


teaches  us  not  to  trust  our  instincts  too  implicitly.  Figure  9.5 
shows  successive  logarithmic  power  spectra  obtained  from  longitudinal 
vibrations  of  a  steel  rod,  approximately  one  meter  long  and  1.5  cm. 
in  diameter  in  its  circular  cross  section.  The  accelerometer  is 
mounted  on  one  end  and  the  other  end  is  struck  sharply  with  a  heavy 
steel  hammer.  Through  the  first  three  modes  the  damping  rates  grow 
quite  rapidly  but  then  settle  down  to  nearly  constant  rates,  exhibit¬ 
ing  no  demonstrable  dependence  on  frequency  from  10  khz.  through  50 
khz.  .  This  experiment  has  been  repeated  many  times  with  much  the 
same  results  in  all  cases. 
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In  reality,  this  result  is  quite  consistent  with  frequency  pro¬ 
portional  damping  in  the  Euler  -  Bernoulli  beam  and  with  our  earlier 
suppositions  as  to  the  source  of  the  damping  action.  If  the  same  as¬ 
sumptions  are  made  in  regard  to  the  internal  structure  of  the  rod  as 
were  made  for  the  beam,  the  resulting  model  for  vibrations  in  the  rod 
takes  the  form  of  the  damped  wave  equation  (constant  coefficient  case 
only  shown  here) 
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z(x,t)  denoting  the  longitudinal  displacement  of  the  rod  element 


whose  equilibrium  center  of  mass  is  at  x  and  jD  (t\)  an  interraction 
kernel  of  the  same  sort  as  h(q)  used  earlier  for  the  beam.  Analy¬ 
sis  similar  to  that  used  for  the  infinite  beam  in  Section  4  shows 
that  a  damping  spectrum  of  the  type  observed  in  Fig.  9.5,  i.e.,  quad 
ratio  damping  exponent  versus  frequency  dependence  for  small  o  and 
asymptotically  constant  as  u  -»  <*,  results  for  the  infinite  rod  and, 
correspondingly ,  similar  discrete  behavior  is  to  be  expected  for  a 
finite  rod.  So  we  may  tentatively  accept  (9.1)  as  a  viable  model 
for  longitudinal,  internally  damped  vibrations  of  a  rod. 

With  (9.1)  so  accepted,  we  now  envision  an  idealized  beam  con¬ 
structed,  as  shown  in  Figure  9.6,  from  two  such  rods  rigidly  spaced 
a  distance  2o  apart,  with  a  mass-less,  perfectly  flexible,  but  inex- 
tensible  central  support  structure  forming  the  "back  -  bone"  of  the 
system.  Lateral  displacements  of  the  structure  will  be  described  in 
terms  of  the  graph  of  the  locus  of  the  center  line,  w  =  w(x,t),  as 
shown  in  the  figure. 
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Tip.  9.6 


Using  the  energy  expression 


for  longitudinal  rod  vibrations,  we  find  with  use  of  (9.1)  and  the 
appropriate  boundary  conditions,  e.g., 

-5z(0,t)  =  -^z(L,t)  =  0  (9.2; 


in  the  case  of  free  ends,  that  for  smooth  solutions  we  have 


dE(z,^)  = 
tb.  at 


i(x-C)  ' 


9z  ,v  +.  i 
-  V  s- ,  t } 


^U,t)  d£dx.  (9. 
9t  J 


Let  z(x,t)  denote  the  longitudinal  displacements  of  the  upper 
and  lower  rods  of  the  beam  structure,  respectively.  Then  we  find,  to 
first  or' d e r  i n  ;; ,  t ha t 


2-U,t)  =  - 


L  ,  -SZ  (  X  ,  t  ) 


^2 r  T 

24V,t) 

9x2 


Let  t  ir.g 


, , ,  9v 
(  w , 


E ( z , 

at' 


Vp'!2  dx 


(9.4 


,  with  p  -  2r,  I  =  pp ,  El  =  2;;2E, 


(w,-v;)  = 


1  p  i  -rW  i  +  I  I  - 
.  K^at-;  pl- 


ata>:J 


+  EI 


•**w~  1  2  1  dx 


while  the  dissipation  form  (9.3)  becomes,  with 


h(x-S)  =  2;^(x-U, 


dE(w,^; 

dt  tt 


h  (x-$ 


*■  a  tax 


( x .  t )  - 


-~-U.t)l  d$dx.  (9.6) 
1"  v 


at  ax 


Assuming  I 


very  small  in  comparison  with  both  p  and  EI 


(  i . e .  ,  we  assume  :  very  small  and  E  very  large),  the  condition  that 
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t he  tine  derivative  of  (9.5  )  with  respect  to  t  should  be  (9.6) 
leads  to  the  equation  (3.7),  with  h(x,£)  =  yh(x-$)  there,  and  the 
arc-riated  boundary  conditions.  If  I  is  not  neglected  we  obtain 

the  corresponding  Rayleigh  model  with  damping  term  as  in  the  Euler- 
Bernouili  model  (3.7).  Thus  the  rod  model  (9.1)  is  consistent 
with  our  earlier  developed  beam  models.  Significantly,  as  we  remark¬ 
ed  at  the  end  of  Section  4,  for  the  damped  Rayleigh  model  the  high 
frequency  damping  exponent  predicted  by  that  model  has  the  same  prop¬ 
erties  as  has  been  observed  in  the  laboratory  (Figure  9.5)  and  as  are 
predicted  from,  the  model  equation  (9.1). 
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1 0 .  Concluding  Remarks . 

In  this  paper  we  have  attempted  to: 

(i)  Show  that  laboratory  beams  of  various  types  do  exhibit 

modes  whose  damping  exponents  are  proportional  to  frequency 
at  the  high  end  of  the  spectrum; 

( i i }  Show  that  there  is  a  relatively  simple  mathematical  model 
which  exhibits  comparable  spectral  behavior  and  whose  form 
is  explainable  in  terms  of  plausible  physical  properties  of 
the  beam; 

(iii)  Show  that  numerical  calculations  using  a  spline  based  app¬ 
roximation  of  our  mathematical  model  yield  calculated 
spectral  values  qualitatively  consistent  with  those  ob¬ 
tained  from  laboratory  experimentation  and  predicted  from 
a:.;. lysis  of  the  mathematical  model. 

Some  readers  may  object  that  the  motivation  of  the  damping  term 
in  our  equation  (3.7)  is  based  on  a  heuristic,  plausibility  argument 
rather  than  on  the  accepted  modus  cperardi  of  introducing  an  approp¬ 
riate  constitutive  equation  and  deriving  the  dynamics  from  that  point 
of  departure  f’®"1  About  all  we  can  say  in  defense  of  ourselves  on  this 
point  is  that  we  are  probably  not  the  best  person  equipped  to  do 
this;  we  have  presented  sufficient  motivation  to  allow  others  who 
specialize  in  continuum  mechanics  to  carry  out  such  a  program  if  they 
are  inclined  to  do  so.  Along  the  same  lines,  we  also  realize  that 
cur  interpretation  of  the  action  of  the  damping  term  in  terms  of  fib- 
e:  structure  of  the  beam,  or  comparable  action  of  extended  crystal 
structures,  may  well  turn  out  to  be  inadequate  in  the  light  of  other 
information  about  the  manner  in  which  damping  forces  act  in  elastic 
structures . 


We  do,  without  apology,  insist  that  laboratory  measurement  of 
damping  rates,  though  necessarily  flawed  to  some  extent  by  exogenous 
influences,  which  we  have  taken  some  pains  to  analyze  here,  must  be 
the  final  arbiter  in  deciding  between  putative  damping  mechanisms. 
Even  spotlessly  correct  theoretical  reasoning  cannot  ultimately 
take  precedence  -  though  it  could  indicate  the  need  for  more  careful 
laboratory  work.  K  ;r:dreds  of  experiments  on  our  part  do  indicate 


frequency  proportional  damping  at  the  high  end  of  the  spectrum  and 
some  sort  of  quadratic  behavior,  or  at  least  the  dominance  of  a  dif¬ 
ferent  proportionality  constant,  at  lower  frequencies.  We  invite 
other  interested  researchers  to  carry  out  similar  experimentation 
and  model  development  complementing,  or  perhaps  even  replacing,  that 
reported  here. 

The  reality  of  near  quadratic  dependence  of  the  damping  expon¬ 
ent  on  frequency  at  the  low  end  of  the  spectrum  seems  to  us  to  be  a 
very  important  question.  Just  citing  the  models  introduced  in  this 
article,  whether  or  not  this  phenomenon  is  real  appears  to  be  the 
decisive  factor  in  choosing  the  model  which  we  have  proposed  in  pre¬ 
ference  to  "square  root"  damping  as  in  (2.3)  which  we  have,  for  the 
present  at  least,  rejected  as  not  adequately  explainable  in  physical 
ter::.  -  -  though  we  have  conceded  that  it  could  be  justified  in  terms 
• f  distances  covered  by  travelling  waves  of  a  given  frequency.  If 
it  ecu  Id  be-  conclusively  demonstrated  that  low  frequency  modes  are 
damps.-.:  a*  rates  which  bear  the  same  relationship  to  frequency  as  the 
i  :  •  rates  associated  with  higher  frequency  modes  then  the  whole 

r  ■  '£  j  process  would  have  to  be  reconsidered  in  that  light.  To 
-'l'.-  "'•.I.  mn-ter  it  would  seem  to  be  necessary  to  conduct  experi- 
w :  * :  **e  ry  long  beams  having  a  large  number  of  low  frequency 

I‘  seems  tr  us  that  the  frequency  independent  damping  of 

higher  m:des  apparently  applying  to  longitudinal  rod  vibrations  is 
a  sri * i ru  1  test  of  the  theory.  It  is  particularly  significant  that 
this  mode  of  damping  is  consistent  with  the  very  different  sort  of 
damping  observed  in  the  beam  and  it  is  intriguing  that  the  mode  of 
dairq  inj  which  apparently  applies  for  the  rod  is  what  is  necessary 
t:  preserve  controllability  by  means  of  finite  dimensional  controls. 
If  frequency  proportional  damping  exponents  did  apply  to  longitudinal 
vibrations  of  the  rod,  those  vibrations  would  be  uncontrollable  in 
either  the  sense  standard  for  the  one  dimensional  wave  equation  [  L  ] 
or  that  which  applies  to  parabolic  processes  [  M  ].  The  whole  area  is 
a  most  intriguing  one.  If  this  article  should  serve  to  encourage  mer 
researchers  to  enter  the  area,  even  if  only  to  refute  our  contribu¬ 
tions,  this  article  will  have  served  its  purpose. 
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